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Abstract 

Calculations  of  the  electronic  mobility  and  drift  velocity  were  carried  out  for  bulk 
GaN  and  AlGaN-GaN  heterojunctions  based  on  a  Monte  Carlo  approach.  Very  good 
agreement  with  available  experiments  was  obtained,  and  the  calculations  yielded  a  set  of 
best-fit  transport  parameters.  The  calculations  were  then  extended  for  the  large-signal 
nonlinear  response  characteristics  of  HEMTs  with  particular  emphasis  on 
intermodulation  distortion  (IMD).  The  results  demonstrated  an  optimal  operating  point 
for  low  IMD  at  reasonably  large  output  power  due  to  a  minima  in  the  IMD  curve.  An 
increase  in  dynamic  range  with  temperature  has  been  predicted,  due  to  a  relative 
suppression  of  interface  roughness  scattering.  The  results  indicate  increasing  the  mole 
fraction  for  the  barrier  layer,  reducing  the  transit  length,  and  introducing  a  thin  AIN 
interfacial  layer  for  suppressing  real  space  transfer  performance  enhancements. 

Simulations  were  also  carried  out  to  probe  possible  enhancements  in  Gun  effect 
oscillators.  It  has  been  shown  that  repetitive  structures  with  serial  segments  to  fashion  a 
"multiple  domain"  device  could  lead  to  significant  improvements  in  output  power  over 
conventional,  single-transit  structure,  and  so  such  multiple  GaN  diodes  merit  serious 
experimental  study.  Monte  Carlo  studies  of  GaN  HEMTs  also  probed  the  effects  of 
scattering  due  to  edge  dislocations  strains.  For  self-consistency,  numerical  solution  of 
Schf'odinger,  Poisson,  and  charge  balance  equations  were  used  for  the  eigen-functions. 
Electron  mobilities  predictions  around  1.711  x  10^  cm^A's,  are  in  close  agreement  with 
reported  data.  At  the  highest  dislocation  density  of  10^®  cm"^  ,  mobility  reductions  of 
16.8  percent  and  8.6  percent  are  predicted  for  77  K  and  300  K,  respectively. 
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1.  INTRODUCTION 

Over  the  past  five  years  or  so,  the  nitride  material  system  has  been  the  focus  of 
intense  research.^'^  Advantages  of  these  direct,  large-bandgap  materials  include  the 
ability  to  sustain  large  electric  fields  (~  10®  Vm‘\  which  bodes  well  for  device 
downscaling  and  higher  breakdown  voltage),  high  saturation  velocities  (~  1.5  x  10  ms'  ), 
lower  generation  noise,  radiation  hardness,  and  high-temperature  operation.  GaN-based 
electronic  amplifiers  are  projected  to  have  significant  improvement  over  silicon  devices.^ 
The  nitrides  have  also  shown  promise  as  emitters  and  detectors  '  ,  Bragg  reflectors, 
light-emitting  diodes  (LEDs),^"^  sensors  for  jet  and  automobile  engines,^^  and  space-based 
operation  in  the  solar  blind  region  (~  260  -  290  nm). 

The  high-field  drift  velocity  of  bulk  GaN  is  larger  than  that  of  GaAs^®’^^  due  to 
the  higher  inter-valley  separation  and  the  larger  optical  phonon  energy.  GaN  lends  itself 
to  hetero-structure  fabrication,  and  the  presence  of  a  strong  internal  polarization  can 
create  very  high  sheet  carrier  densities  (  >  10^^  m'^).  Experimental  demonstrations  of 
large  radio-frequency  (RF)  power  densities,  as  measured  in  Watts  per  millimeter  of  the 
gate  periphery,  have  already  been  made^®'^\  The  higher  thermal  conductivity  of  GaN 
relative  to  GaAs  and  use  of  SiC  substrates^^’^®  should  help  alleviate  the  thermal 
management  issue.  These  aspects  collectively  enhance  the  prospects  of  GaN  for 
microwave  power  amplification,  particularly  at  the  Z-band  and  higher  frequencies.^'^'^^ 

Regarding  electron  transport  properties,  the  5-valley  effective  mass  of  GaN  is 
0.22,  which  is  about  three  times  higher  than  that  of  GaAs.  Hence,  the  mobility  in  bulk 
GaN  is  expected  to  be  inferior  to  GaAs,  and  this  has  been  shown  in  experimental  work. 
However,  the  high-field  drift  velocity  of  bulk  GaN  is  larger  than  that  of  GaAs.^’*^  This 
is  due  to  the  higher  inter-valley  separation  and  the  larger  energy  of  optical  phonons.  The 
former  helps  keep  a  higher  fraction  of  electrons  within  the  lower  mass  valley,  while  the 
latter  reduces  the  net  role  of  phonon  scattering  by  increasing  the  emission  threshold.  A 
more  striking  advantage  for  GaN,  however,  results  upon  carrier  confinement  that  leads  to 
the  creation  of  a  two-dimensional  electron  gas  (2DEG).  Carrier  confinement  is  typically 
engineered  through  the  use  of  single  or  double  modulation-doped  heterostructures^^. 
This,  as  proposed  originally  for  the  GaAs  system,  leads  to  the  physical  separation  of  the 
mobile  electrons  from  their  parent  donors,  thereby  reducing  impurity  scattering.  In  GaN, 
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however,  the  presence  of  a  strong  internal  polarization  arising  from  the  spontaneous  and 
piezoelectric  effects  works  to  create  stronger  band  bending  and  contributes  to  very  high 
sheet  carrier  densities,  even  in  the  absence  of  intentional  doping.  Characterization  of  the 
heterojunction  field-effect  transistors  (HFETs)  fabricated  has  revealed  strong 
enhancement  in  carrier  densities,^  which  should  lead  to  higher  currents.  The  ac  gain  is 
therefore  expected  to  be  large,  and  so  GaN  HFETs  offer  bright  prospects  for  use  in 
microwave  amplification,  particularly  at  the  Z-band  and  higher  frequencies.^^’^® 

Polarization  arises  in  GaN  due  to  its  non  centro-symmetric  atomic  arrangement. 
The  nitrides  exhibit  large  piezoelectric  effects  under  stress  along  the  c-direction  due  to 
the  large  lattice  mismatch.  Piezoelectric  coefficients  are  almost  an  order  of  magnitude 
larger  than  the  traditional  III-V  materials.^^  Such  polarization  can  cause  a  sizeable  Stark 
shift  in  quantum  wells^^’^^  and  works  to  alter  the  interfacial  concentration  of  free  carriers 
in  strained  zincblende  GaN  grown  along  the  (111)  orientation^'^  or  in  pseudomorphic 
wurtzite  material  grown  in  the  (0001)  orientation.^^  Carrier  densities  as  high  as  4  x  10^^ 
cm'^  have  been  reported,  even  in  the  absence  of  doping  within  the  barrier  layer.^®  In 
addition  to  strain-related  electric  fields,  a  large  spontaneous  polarization  exists,  even  at 
zero  strain,  as  pointed  out  by  Bernardini  et  al.^’  Aluminum  nitride,  for  example, 
possesses  a  spontaneous  polarization  that  is  only  about  3  times  smaller  than  that  of  a 
typical  ferroelectric  perovskite.^^  Hence,  interface  fields  as  high  as  2  MV/cm  can  result. 

A  number  of  issues  somewhat  unique  to  the  GaN-based  heterostructure  system 
arise  and  can  be  expected  to  have  important  bearing  on  the  transport  behavior.  The  first 
issue  is  related  to  the  large  polarization.  For  instance,  polarization-related  modifications 
to  the  internal  fields  could  affect  ionization  rates  in  quantum-well  devices^^  or  alter 
electronic  injection  at  the  gate  contacts  of  HFETs.'**^  Furthermore,  carrier  mobility  should 
also  be  affected  by  the  polarization-induced  field  enhancements,  since  carrier  wave 
functions  and  the  scattering  rates  are  directly  modified  in  the  process.  In  addition, 
interface  roughness  scattering  is  expected  to  increase  due  to  the  close  proximity  of  the 
2DEG,  with  the  barrier  layer  at  the  higher  electric  fields.  The  interface  roughness 
scattering  is  typically  quantified  on  the  basis  of  weak-perturbation  theory,  in  terms  of  two 
parameters:  the  root-mean-square  value  of  fluctuations  at  the  interface  ),  and  the 
correlation  length  L  between  fluctuations.  These  parameter  values  are  not  well  known  for 
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AlGaN/GaN  HFETs  and  need  to  be  deteraiined.  It  is  also  unclear  whether  weak 
perturbation  theory  would  be  adequate  for  accurately  predicting  the  transport  behavior. 
Screening  effects  arising  from  the  higher  electron  density  also  need  to  be  considered  in 
evaluating  the  transport  properties.  These  could  exhibit  spatial  inhomogeneity  across  the 
interface.  Next,  the  degeneracy  effects  associated  with  the  higher  electron  density  need 
to  be  carefully  evaluated  because  of  the  high  2DEG  density.  There  have  been  no  studies 
into  the  role  and  extent  of  Pauli  exclusion  on  the  mobility  in  GaN  HFETs,  to  the  best  of 
our  knowledge.  A  third  point  relates  to  the  impact  of  material  parameters  and  geometry 
on  the  electrical  response.  For  example,  changes  in  the  mole  fraction  of  the  AlGaN 
barrier  layer,  the  spacer  thickness,  or  the  barrier  width  all  alter  the  polarization  field,  the 
carrier  density,  and,  thus,  the  electrical  response  characteristics.  It  is  important  to  obtain 
quantitative  predictions  in  order  to  optimize  the  device  operating  characteristics.  Finally, 
numerical  evaluations  of  the  drift  velocity  as  a  function  of  the  applied  gate  voltage  and 
longitudinal  field  need  to  be  carried  out  to  ascertain  the  nonlinearity  in  response  and  the 
high-temperature  behavior.  Changes  in  the  interface  field  with  gate  bias,  for  example,  are 
expected  to  produce  nonlinear  variations  in  drift  velocity.  Such  nonlinearity  can  lead  to 
possible  harmonic  generation,  intermodulation  and  mixing  of  time-dependent  signals 
applied  to  the  gate  of  HFET  structures. 

One  of  the  most  important  requirements  for  microwave  power  amplifiers  is  the 
level  of  nonlinear  distortion  and  their  inter-modulation  (IM)  behavior.  Nonlinearity  in  the 
output  current  response  with  applied  gate  voltage  results  from  the  combined  nonlinear 
behaviors  of  carrier  drift  velocity  and  the  channel  density.  The  scattering  rates,  which 
control  drift  velocity,  are  nonlinear  functions  of  energy.  For  example,  interface 
roughness  scattering  exhibits  a  nonlinear  increase  with  gate  bias"^^  due  to  the  increased 
proximity  of  the  two-dimensional  electron  gas  (2DEG)  with  the  barrier  layer.  The  gate 
voltage  (Vg)  dependent  variations  in  the  2DEG  density,  also  affect  scattering  through 
changes  in  the  screening  of  polar  interactions  and  the  degeneracy  level  Carrier  density 
dependence  on  the  mole  fraction,  barrier  thickness,'^^  and  real  space  transfer  into  AlGaN, 
which  depends  on  both  the  gate  and  drain  bias  levels,  all  work  to  modify  the  transport 
characteristics  further. 
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Since  the  satellite  valley  effective  electronic  masses  of  both  zincblende  and 
wurtzite-phase  GaN  are  higher  than  in  the  lowest  conduction  bands,  the  material  appears 
to  show  promise  as  a  large-signal  oscillator  based  on  the  transferred  electron  effect.'^^ 
The  presence  of  an  inflection  point  in  the  3i-valley  bandstructure  of  the  zincblende  phase 
has  also  been  suggested  as  a  source  of  negative  differential  resistance  (NDR)'^'^  .  As 
compared  to  the  traditional  GaAs  and  InP-based  negative  differential  resistance  (NDR) 
devices, the  conversion  efficiency  and  power  handling  capacity  should  be  superior. 
Despite  the  potential,  there  have  not  been  many  studies  on  this  subject  to  the  best  of  our 
knowledge.  The  only  reports  appear  to  be  by  Alekseev  and  Pavlidis"^^  for  wurtzite  GaN 
based  on  transient  hydrodynamic  simulations,  and  by  Zhao  et  al.^^  for  the  zincblende 
phase.  A  full  Monte  Carlo  treatment  of  the  time-dependent  transferred  electron  problem, 
is  more  accurate  and  desirable  than  the  hydrodynamic  scheme  Previous  studies  have 
often  used  the  steady-state  velocity-field  characteristics'^  and  relied  on  somewhat  older 
band  structure  data.  For  example,  previous  reports  were  based  on  3i-32  and  M-L 
intervalley  separations  of  1.9  eV  and  2.1  eV,  respective! More  recent  calculations^^ 
suggest  a  3] ,  M-L,  33 ,  Li ,  Ai  ordering,  with  the  two  lowest  valleys  2.2717  eV  and  2.4 
eV  above  3i  .  In  view  of  this,  it  is  perhaps  important  to  re-examine  the  large-signal 
response  characteristics  of  GaN-based  oscillators  based  on  recent  data. 

2.  MONTE  CARLO  SIMULATIONS 

2.1  Bulk  GaN  Material 

The  simulation  scheme  for  bulk  wurtzite  GeiN  is  discussed  first.  Details  of  the 
Monte  Carlo  approach  are  similar  to  those  described  previously.'^  A  three-valley, 
nonparabolic  conduction  band  model  was  used  for  the  transport  calculations  in  bulk  GaN. 
For  wurtzite-phase  GaN,  the  minimum  is  located  at  the  >point.  The  satellite  valleys  are 
at  the  f/-point  that  is  two-thirds  of  the  way  between  the  L-  and  M-  symmetry  points.  The 
next  highest  valley  is  located  at  the  33-point.  Thus,  two  ^  -valleys  and  the  six  equivalent 
[/-valleys  were  considered.  A  three-dimensional  view  of  the  Brillouin  zone  including 
points  within  the  irreducible  wedge  is  shown  in  Fig.  1.  The  effective  masses  have  been 
calculated  to  be  0.2283  and  0.1846  along  the  perpendicular  and  parallel  directions, 
respectively,  with  a  non-parabolicity  coefficient  of  about  0.37  eV'^  The  satellite  valleys 
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Fig.  1.  (a)  Brillouin  zone  schematic  including  the  various  symmetry  points  within  the 
irreducible  wedge,  (b)  View  of  the  full  brillouin  zone  including  the  irreducible  wedge. 

are  at  the  l/-point  that  is  two-thirds  of  the  way  between  the  L-  and  M-  symmetry  points, 
with  an  energy  separation  of  2.2717  eV.  According  to  the  calculations  by  Goano  et  al.^\ 
the  next  highest  valley  is  located  at  the  33-point,  2.4  eV  above  33 .  Finally,  the  valleys  at 
the  Li  and  Ai  points  are  located  2.5  eV  and  2.6  eV,  respectively  above  the  conduction 
band  minima.  The  effective  masses  for  the  M-L  minima  are  anisotropic.  As  given  in 
Table  I  the  mass  is  3.0375  (along  3-A1) ,  0.3158  (along  K-H),  and  0.3858  (along  the  M 
and  L  points).  The  corresponding  nonparabolicity  coefficients  also  depend  on  the 
direction,  and  are  reported  to  be:  -1.0325  eV^  (along  3-Ai)  ,  0.2164  eV‘^  (along  K-H),  - 
0.174  eV^  (along  M)  and  0.9686  eV^  (along  L).  Such  anisotropic  features  were  included 
in  the  present  Monte  Carlo  approach.  For  the  33  valley,  the  effective  masses  along  the 
parallel  and  perpendicular  directions  were  taken  to  be  3.6227  and  0.2856,  respectively. 
The  nonparabolicity  factors  are  given  to  be  4.4554  parallel  to  the  c-axis  and  0.22  along 
the  perpendicular  directions. 

The  scattering  mechanisms  considered  were  ionized  impurity  based  on  the 
Brooks-Herring  approach,  acoustic  deformation  potential  scattering,  polar-optical 
interactions,  and  intervalley  deformation  potential  processes.  The  piezoelectric  scattering 
was  excluded,  since  it  has  been  shown  to  be  negligible  at  temperatures  of  300  K  and 
beyond.  Nonequivalent  intervalley  scattering  events  were  taken  into  account  amongst  all 
three  valley  types,  governed  by  a  single  deformation  potential  and  phonon  energy.  A  65- 
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Table  I  Parameters  used  in  the  Monte  Carlo  calculations  for  GaN 


PARAMETER 

VALUE 

Lattice  constants  in  meter  (a  ,  c) 

3.189  X  10'^“ ,  5.185  X  10'^*^ 

Dielectric  constants  k(0) ,  k(4) 

9.5 ,  5.35 

Density  (gm  cm'^) 

6.095 

Acoustic  velocity  (ms  ■^) 

4.33  X  10^ 

Effective  mass  {Ti,U-  and  T3 -valley) 

0.21 ,  0.25  ,  0.4 

Valley  separation  (eV) 

0.0,  1.95,2.1 

Nonparabolicity  factors  for  the  three  valleys  (eV^) 

0.19 , 0.1 ,  0.0 

Longitudinal  optical  phonon  energy  (eV) 

0.092 

Intervalley  phonon  energy  (eV) 

0.065 

Acoustic  deformation  potential  (eV) 

8.0 

Zero-order  equivalent  intervalley  deformation  potential  (eV/cm)  0.5  x  10^ 

Zero-order  nonequivalent  intervalley  deformation  potential  (eV/cm)  1.0  x  10^ 

First-order  intervalley  potential  (eV) 

5.0 

Interface  roughness  correlation  length  (nm) 

1.5 

RMS  fluctuations  of  interface  roughness  (nm) 

0.65 

Aluminum  alloy  fraction 

X 

Spontaneous  polarization  (C/m^) 

-0.052*x  -0.029 

Elastic  constants  C13  and  C33  (GPa) 

5x  +  103  and  -32x  +  405 

AlGaN  alloy  scattering  parameter  (eV) 

0.5 
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meV  phonon  was  assumed  to  adequately  represent  both  nonequivalent  and  equivalent 
intervalley  scattering.  A  92-meV  phonon  value  was  used  for  the  intravalley  longitudinal 
polar  optical  phonon  scattering.  Screening  was  incorporated  with  the  random  phase 
approximation  in  the  long  wavelength  limit.  The  value  of  the  inverse  screening  length  3 
was  obtained  as:  .^t)  =  [n  ^  l,kBTe(fyf  '^  where  n  is  the  electron  density,  e  the  electronic 
charge,  ks  the  Boltzmann  constant,  and  Jg  the  effective  electron  temperature.  This 
effective  temperature  was  evaluated  at  each  time  step  of  the  Monte  Carlo  simulation  with 
the  following  equation: 

Te(t)  =  [2/(3kB)]  3  ,  //  {  0.5  mi  <v^>/  -  0.5  mi  }  •  (1) 

In  the  above,  //,  mi,  <v^>/  and  v^di  refer  to  the  electron  fraction,  effective  mass,  mean- 
square  velocity,  and  drift  velocity,  respectively,  for  the  band  and  the  index  I  runs  over 
the  three  bands.  The  Tg  represents  an  ad-hoc  parameter  as  used  by  Bhapkar  and  Shur^^ 
and  is  not  a  real  electron  temperature.  Pauli  exclusion  was  included  based  on  a  rejection 
technique  implementation  first  proposed  by  Lugli  and  Ferry.^^  The  material  parameters 
required  for  the  bulk  GaN  simulations  were  generally  taken  from  the  published 
literature,^^’^'*  with  minor  adjustments. 

2.2  GaN-AlGaNHEMTs 

Simulations  of  parallel  transport  for  the  2DEG  in  GaN-AlGaN  HFET  structures 
were  also  carried  out  based  on  the  kinetic  Monte  Carlo  approach.  As  is  well  known, 
mobile  carriers  are  confined  within  a  shallow  inversion  layer  at  the  interface. 
Confinement  then  leads  to  the  following  consequences:  (i)  electron  quantization  with 
subband  formation,  leading  to  a  two-dimensional  flow;  (ii)  changes  in  the  transverse 
wavefunctions  that  slightly  lower  the  electron-phonon  scattering  rates,  but  increase  the 
number  of  inter-  and  intra-subband  processes;  (iii)  a  relative  decrease  in  the  influence  of 
remote  ion,  oxide  charge,  and  depletion  dopant  ion  scattering  due  to  the  increased 
physical  separation;  and  (iv)  the  inclusion  of  an  additional  interface  roughness  scattering, 
which  imparts  a  strong  transverse  electric  field  dependence  to  the  channel  mobility,  as 
first  shown  by  Sabnis  and  Clemens^^  for  silicon.  Physically,  higher  transverse  fields 
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enhance  the  electronic  confinement  by  drawing  them  closer  to  the  interface,  which  leads 
to  greater  interface  roughness  scattering  and  lower  mobility.  The  details,  however,  are 
complicated,  since  the  relative  subband  occupancy  and  spatial  electron  distribution  within 
each  subband  change  with  the  transverse  field  by  different  amounts.  Furthermore,  as 
channel  densities  are  enhanced  with  transverse  field,  free-carrier  screening  increases  and 
should  work  to  lower  the  effects  of  polar  and  Coulombic  interactions,  thereby  moderating 
the  mobility  reductions. 

For  simplicity,  it  has  been  assumed  here  that  the  electron  density  within  the 
conduction  channel  is  uniform.  Extensions  to  account  for  possible  nonuniformities  can 
be  incorporated  by  using  the  procedure  outlined  by  Ravaioli  and  Ferry.^®  Details  of  our 
inversion-layer  Monte  Carlo  implementation  and  procedure  can  be  found  elsewhere. 
Electron  confinement  was  treated  on  the  basis  of  a  simple  triangular  well 
approximation.^®’^^  In  keeping  with  the  recent  experimental  work  on  GaN 
heterostructures,  the  c-axis  was  taken  to  be  perpendicular  to  the  GaN-AlGaN  interface 
(and  hence,  parallel  to  the  growth  direction),  resulting  in  a  single  transverse  effective 
mass  and  one  set  of  subbands.  Fang  and  Howard^”  variational  wavefunctions  were  used 
for  computations  of  the  scattering  rates.  Although  more  sophisticated  wavefunctions 
derived  from  a  solution  of  the  Poisson  and  Schrodinger  equations  could  have  been  used, 
the  present  approach  was  chosen,  given  the  existing  uncertainty  in  the  material 
parameters  and  the  need  to  implement  a  relatively  fast  computational  scheme.  Besides, 
the  variational  method  is  known  to  yield  reasonably  accurate  results  and  to  exhibit  the 
correct  trends.  A  two-subband  Monte  Carlo  model  was  used,  that  included  electron 
interactions  with  acoustic  modes  via  the  deformation  potential,  polar  optical  phonon 
interactions,  zero-  and  first-order  intervalley  deformation  potential  scattering,^^  interface 
roughness,^^  and  interface  ion  interactions.^^  Higher  subbands  were  ignored  in  the  present 
context  of  low  fields,  as  were  considerations  of  interface  phonon  modes.  Similarly,  real- 
space  transfer  and  electron  flow  in  the  AlGaN  barrier  layer  was  also  neglected.  However, 
the  role  of  AlGaN  on  the  transport  was  indirectly  taken  into  account  by  including  an  alloy 
disorder  scattering  mechanism.  The  rationale  is  similar  to  that  used  by  Hsu  and 
Walukiewicz.  Screening  of  the  polar  optical  and  interface  ion  interactions  was 
evaluated  within  the  random  phase,  zero-frequency  approximation^'^  that  requires  the 
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evaluation  of  a  wavevector-dependent  form  factor.  Finally,  both  intra-  and  interband 
transitions  were  accounted  for  in  the  numerical  transport  simulations.  Treatment  of  the 
interface  roughness  scattering  was  based  on  a  simple  Gaussian  model  for  the 
autocovariance  function,^®  rather  than  the  exponential  dependence  used  by  Goodnick  et 
al.  for  silicon.^^  More  complicated  treatments  of  discrete  AlGaN  islands  at  the  interface 
have  been  suggested  to  include  localization  effects, but  were  not  used  here  for 
computational  simplicity.  The  interface  roughness  scattering  is  most  simply  treated  in  the 
Born  approximation  by  assumption  that  the  roughness  can  be  regarded  as  a  weak 
perturbation.  This  is  then  characterized  in  terms  of  two  parameters,  namely,  the 
correlation  length  L  and  root-mean-square  (rms)  height  )  of  the  roughness  fluctuations. 
This  model  has  previously  been  shown  to  be  adequate,  yielding  L  and  )  values  of  about 
2.2  nm  and  0.2  nm  for  silicon®^  and  2.3  nm  and  0.8  nm  for  SiC.^^  The  corresponding 
parameters  for  the  GaN-AlGaN  system,  however,  are  not  known,  and  have  been  obtained 
here  through  a  best-fit  procedure.  Finally,  degeneracy  is  expected  to  play  a  role  in  the 
GaN  2D  system,  because  of  the  large  carrier  densities  that  have  been  observed 
experimentally.  This  effect,  which  essentially  reduces  the  scattering  rate  due  to  a  finite 
occupancy  probability  of  the  final  state,  was  taken  into  consideration. 

2.3  Polarization  Effects  in  HEMTs 

As  alluded  to  earlier,  the  nitrides  lack  inversion  symmetry  and  exhibit  strong 
piezoelectric  effects.  Such  noncentro-symmetric  compounds  exhibit  two  different 
sequences  of  their  atomic  layering.  For  GaN  heterostructures  with  the  growth  direction 
normal  to  the  {0001}  basal  plane,  the  atoms  are  arranged  in  bilayers.  These  bilayers 
consist  of  two  closely  spaced  hexagonal  layers,  one  formed  by  cations  and  the  other  by 
anions,  leading  to  polar  faces.  Thus,  a  basal  surface  could  be  either  Ga-  (or  N-)  faced, 
corresponding  to  either  Ga  (or  N)  atoms  on  the  top  position  of  the  {0001}  bilayer.  The 
two  (0001)  and  (0001)  surfaces  are  not  equivalent,  and  lead  to  different  net  polarizations 
and  internal  charge  densities.  For  clarity,  an  HFET  with  its  equilibrium  energy-band 
diagram  is  shown  in  Figs.  2(a)-2(c)  for  the  Ga-  and  N-polarity.  The  GaN  is  relaxed, 
while  the  AlGaN  has  been  assumed  to  be  under  tensile  strain.  The  spontaneous 
polarization  for  AlGaN  is  larger  in  magnitude  than  that  for  GaN.  The  spontaneous  and 
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Fig.  2.  Equilibrium  band  diagrams  of  a  GaN-AlGaN  HFET:  (a)  undoped  GaN  with  Ga- 
polarity,  (b)  n-type  GaN  with  N-polarity,  and  (c)  p-type  GaN  with  N-polarity. 
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piezoelectric  polarizations  support  each  other,  and  a  large  interface  charge  is  created.  In 
Fig.  2a,  the  resulting  potential  slopes  toward  the  interface  from  both  sides,  and  helps 
drive  mobile  electrons  toward  the  channel.  There  is  thus  no  specific  need  to  dope  the 
AlGaN  or  have  a  spacer  layer.  The  metallic  reservoir  and  carrier  generation  processes 
can  furnish  the  requisite  electrons  into  the  GaN.  The  situation  for  a  N-polarity,  on  the 
other  hand,  is  different,  as  shown  in  Figs.  2b  and  2c.  Both  n-  and  jn-type  doping  have 
been  considered,  but  the  net  result  is  the  same.  The  energy  diagram  predicts  an 
accumulation  of  holes  at  the  interface  rather  than  electrons.  In  any  case,  it  becomes  clear 
that  details  of  the  surface  polarity  need  to  be  taken  into  account,  and  the  most  effective 
scenario  for  a  GaN  HFET  appears  to  be  the  Ga-polarity  that  would  create  a  large  electron 
pool  without  having  to  incorporate  any  dopant  impurities. 

The  piezoelectric  polarization  Ppe  can  be  calculated  from  the  piezoelectric 
constants  633  and  e3i  as:  Ppe  =  633  %  +  €31  %)  ,  where  %  is  the  strain  along  the  c- 

axis  and  the  in-plane  strains  '"fx  and  'fy  are  equal  and  assumed  to  be  isotropic.  The  strains 

can  be  expressed  in  terms  of  changes  in  the  lattice  constant  parallel  («„)  and 
perpendicular  (Cg)  to  the  c-axis  as: 

~  /  Og  ,  (2U) 

%  =  [c(x)  -Cg]l  Cg  =  -2  [Ci3(x)l  C33(x)]  {[a(x)  -Ug]/  ag}  .  (2b) 

In  the  above,  and  C33(x)  are  the  elastic  constants,  x  is  the  aluminum  mole  fraction 

in  the  AlxGai-xN  material,  and  a(x)  is  the  in-plane  lattice  constant.  Values  of  the 
piezoelectric  coefficients  e33(x)  and  as  well  as  the  elastic  constants  Ci3(x)  and 

C33(x)  for  the  GaN-AlxGai-xN  system  can  be  expressed  in  terms  of  the  mole  fraction  x 
as:'^ 

e33(x)  =  (0.73  x  + 0.73)  Cm‘^  and  633  (x)  =  (-0.11  x  -  0.49  )  Cm'^  ,  (3a) 

C33  (x)  =  (-32  X  +  405  )  GPa  and  C13  (x)  =  (5  x  +  103  )  GPa  .  (3b) 

Finally,  the  spontaneous  polarization  Psp(x)  for  this  material  system  is  expressed  as: 
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Psp(x)  =  (-  0.052  a:  -  0.029 )  Cm'^  . 


(4) 


Hence,  the  magnitude  of  the  polarization-induced  interface  sheet  charge  density  ®s(x)  in 
Cm'^,  for  an  undoped,  pseudomorphic,  AlGaN-GaN  single  heterostructure  turns  out  to  be 

^x)  =  0.0483x{[(5x+103)/(-32x  +  405)]  (0.73x+0.73)  +  (0.1b:+0.49)  }  +  0.052  a:  .  (5) 

This  interface  charge  works  to  alter  the  interface  electric  field,  and  hence,  affects  the 
channel  density,  the  interface  roughness  scattering,  the  free  carrier  screening,  and  the 
channel  degeneracy.  The  channel  electron  charge  density  qns(x),  is  then : 

q«i(x)  =  ^s(x)  -  [,o,(-^)/d]  {Nb(x)  +  Ef(x)  -  )Ec(x)  }  ,  (6) 

with  ,0  being  the  permittivity  of  free  space,  ,(x)  the  dielectric  constant,  d  the  AlxGai.xN 
barrier  thickness,  Nb(x)  the  Schottky-barrier  height,  Ef(x)  the  Fermi  level  with  respect  to 
the  GaN  conduction  band-edge,  and  )Ec(x)  the  conduction  band  offset.  Using  values 
reported  in  the  literature,^®’^^  leads  to  : 


,(x)  =  -  0.5  X  -H  9.5  ,  (7a) 

Nb(x)  =  (1.3  x  + 0.84)  ,  (7b) 

)Ec(x)  =  0.7[Eg(x)-Eg(0)]  ,  and  (7c) 

Ef(x)  =  {[9BSq^  tisix)]  /[8,o  ,(x)  {8m*(x)}°-^]  +  [BE^  ni(x)]/m*(x)  .  (7d) 

The  effective  mass  m*(x)  ~  0.22  x  9.1  x  10'^^  Kg,^°  while  the  energy  gap  Eg(x)  is 
measured  in  eV  to  be^^ : 


Eg(x)  =  6.13  X  +  3.42  (1-x)  -  x  (1-x)  .  (7e). 

3.  GaN  HEMTs  -  ROLE  OF  DISLOCATIONS 

Epitaxial  layers  of  GaN  can  contain  a  high  density  of  dislocations,  which  form 
because  of  the  absence  of  a  suitably  lattice-matched  substrate  for  growth.  Both  sapphire 
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and  silicon  carbide  substrates  have  such  problems.  These  dislocations  are  predominantly 
oriented  parallel  to  the  c  axis  of  the  material. Dislocations  densities  in  GaN  layers 
grown  on  sapphire  or  SiC  substrates  range  between  10^-10^  cm’^  One  of  the  important 
effects  of  dislocations  that  have  been  reported  in  the  GaN  arena  is  on  carrier  scattering^"'' 

'jc 

due  to  the  creation  of  localized  charge  centers.  For  example,  room  temperature 
electron  mobilities  were  shown  to  decrease  from  about  600  cm  fWs  to  150  cm  A^s  as  the 
dislocation  density  increased  from  4  x  10^  to  2  x  lO'®  cm'^  in  bulk  GaN  films.^®  Using  a 
Boltzmann  transport  model  with  dislocation  scattering  treated  in  terms  of  a  screened 
interaction  between  mobile  electrons  and  line-charges,  Look  and  Sizelove’^  were  able  to 
fit  the  bulk  GaN  mobility  data  rather  well.  The  electron  mobility  vs  density  was  shown 
to  have  a  local  maxima  and  decrease  with  density  at  low  concentrations,  due  to 
reductions  in  the  screening  of  dislocation  scattering.  Another  unusual  transport 
phenomena  in  bulk  GaN,  attributed  to  dislocations,  was  that  of  a  mobility  collapse  at  low 
carrier  densities.^^  In  a  theoretical  analysis,^^  it  was  shown  that  this  sharp  density 
dependent  transition  between  regions  of  low-  and  high-mobilities  could  arise  from  a 
collective  effect  of  dislocation  walls  leading  to  the  formation  of  internal  potential 
barriers.  Thus,  the  role  of  dislocations  on  the  transport  in  bulk  GaN  appears  to  be  well 
recognized  and  shown  to  have  a  strong  impact  on  carrier  mobilities. 

By  contrast,  the  role  of  dislocations  on  the  transport  of  a  two-dimensional  electron 
gas  (2DEG)  has  not  been  studied  in  detail.  Some  reports  of  carrier  mobilities  in  high 
electron  mobility  transistors  (HEMTs)  indicate  that  dislocation  scattering  may  not  play 

Qr\ 

such  a  significant  role.  This  suggests  that  the  role  of  screening,  which  is  known  to  be 
somewhat  different  in  two-dimensional  systems,  might  not  be  quite  as  strong.  Part  of  the 
resolution  may  also  lie  in  the  details  of  the  mobility  data.  For  instance,  the  mobility 
values  at  300  Kelvin  were  shown  to  decrease  only  mildly,  with  variations  in  the  2DEG 
concentration.  The  mobilty  was  reported  to  vary  from  about  1550  cmA^s  to  1000 
cm^A^s.*®  The  range  of  densities  for  these  cases,  however,  was  roughly  around  lO'^  cm'^, 
and  hence,  fairly  large  on  an  absolute  scale.  By  comparison,  channel  densities  in  the 
well-known  GaAs-heterostmcture  systems,  are  seldom  as  high  as  lO'^  cm'^  .  Screening, 
which  might  be  responsible  for  such  mobility  effects,  is  quantitatively  quite  complicated 
in  a  multi-subband  2D  system.  Its  evaluation,  even  in  the  random  phase  approximation 
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(RPA),  requires  an  inversion  of  a  dielectric  matrix.^^  For  the  simplest  situation  of  single 
subband  static  screening  under  the  RPA  and  in  the  long-wavelength  limit,  the  dielectric 
function  can  be  characterized  by  a  wavevector  independent  inverse  screening  length  3. 
This  inverse  screening  length  is  independent  of  the  carrier  density  at  low  temperatures*^, 
and  proportional  to  the  2DEG  concentration  under  non-degenerate  conditions.*^ 
Expressions  for  3  under  these  two  limiting  conditions  are: 

3  (low  temperature)  ~  {  m*  e^  }/  [BZ^  (,i  +  ,2 )  ]  ,  (8a) 

and  ,  3  (non-degenerate)  ~  [  n2D  ^  ]/  [(,1  +  ,2  )  keT  ]  ,  (8b) 

where,  ,1  and  ,2  are  the  AlGaN  and  GaN  material  permittivities,  m*  the  effective 
electronic  mass,  n2D  the  2DEG  density,  ks  the  Boltzmann  constant,  and  T  the  temperature 
in  Kelvin. 

Clearly,  at  low  electron  concentrations,  and/or  high  temperature  operation  as  would 
practically  occur  for  GaN  HEMTs,  eqn.  (8b)  would  be  more  appropriate.  In  the  GaN- 
based  experimental  devices,  however,  the  2DEG  density  was  relatively  high.  Hence,  the 
screening  length  would  only  have  been  weakly  dependent  on  the  2DEG  density.  This 
might  have  been  a  factor  in  the  apparent  lack  of  a  strong  density  dependent  affect  on  the 
mobility.  Besides,  for  the  room  temperature  experiments,  phonon-scattering  may  be  the 
dominant  interaction,  and  would  then  overwhelm  the  role  of  other  mechanisms  such  as 
dislocation  scattering.  A  more  appropriate  situation  for  probing  the  impact  of 
dislocations  on  2DEG  carrier  transport  would  be  at  low  densities.  Moderate-to-high 
operating  temperatures  might  not  be  constructive,  as  the  phonon  scattering  would 
increasingly  dominate.  The  inverse  screening  at  low  carrier  densities  would  be  reduced, 
and  collective  effects  brought  about  by  dislocations  (such  as  localized  strain  and  charges 
at  the  hetero-interface)  could  prevail.  Such  a  low-density  scenario  is  not  completely 
hypothetical,  since  low  concentrations  from  a  practical  standpoint,  would  arise  for  thinner 
AlGaN  barrier  layers,*'^  devices  with  lower  A1  mole  fractions,*"^  or  during  the  negative 
polarity  excursions  of  the  gate  voltage  for  ac  operation. 

Traditional  evaluations  of  electronic  transport  in  the  presence  of  defects,  have 
considered  Coulomb  interactions  between  mobile  carriers  and  dislocation  charges.*^  The 
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strain  fields  associated  with  defects  have  generally  been  neglected.  However,  recent 
calculations  show  that  no  internal  electric  fields  are  generated  for  screw  dislocations,  and 
that  for  edge  dislocations,  polarization  charge  density  and  electric  fields  are  generated 
only  in  the  near  interfaces'^  This  makes  the  dislocation  scattering  especially  important 
and  germane  to  interface  transport,  as  in  the  case  of  a  two-dimensional  electron  gas 
(2DEG)  in  HEMTs.  In  a  very  recent  report,  the  effect  of  such  edge-dislocation  induced 
strain  fields  on  2DEG  scattering  was  analyzed  to  obtain  an  expression  for  the  strain- 
limited  mobility.  A  useful  closed-form  expression  for  dislocation  mobility  was 
presented,  though  the  results  relied  on  the  relaxation  time  approximation,  and  collective 
contributions  of  other  scattering  mechanisms  (e.g.  carrier-phonon)  were  not  taken  into 
accounted.  In  any  case,  it  was  shown  that  for  a  given  dislocation  density,  the  mobility 
reduction  would  be  the  largest  at  the  smallest  carrier  density  n2D. 

As  is  rapidly  becoming  clear,  the  2DEG  density  n2D  is  controlled  by  several 
factors  such  as  the  thickness  of  the  AlGaN  barrier  layer,  the  presence  of  surface  states 
and  role  of  possible  capping  layers.  *  Hence,  an  explicit  connection  with  the  structure 
and/or  geometry  of  the  HEMT  device  exists.  Fig.  3  shows  an  energy  band  schematic 
for  a  GaN-AlGaN  HEMT  with  a  GaN  cap  layer.  In  general,  the  capping  layer  shown 
may  be  of  another  material  (e.g.  silicon  nitride)  or  might  even  be  totally  absent. 
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Fig.  3.  Energy  band  schematic  for  a  GaN  HEMT  with  an  AlGaN  barrier  layer  and  a  GaN 
capping  layer.  The  surface  potential  is  assumed  fixed  with  negligible  doping  in  AlGaN. 
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Regardless  of  the  details,  the  2DEG  density  depends  on  the  AlGaN  barrier  doping  and  the 
position  in  energy  of  the  interface  donor-like  states  above  the  Fermi  level.  The  latter  is 
controlled  by  the  barrier  and  capping  layer  thicknesses.®"^  It  is,  therefore,  possible  to  have 
situations  (e.g.  HEMTs  with  barriers  having  no  intentional  doping,  but  with  a  capping 
layer)  where  nzo  might  be  low,  with  holes  present  at  the  AlGaN-GaN  cap  interface.  It 
would  also  be  possible  to  lower  the  inversion  charge  upon  the  application  of  a  positive 
gate  potential,  with  the  holes  moving  out  towards  the  2DEG  at  the  GaN-AlGaN  interface. 
This  has  obvious  implications  for  large  signal  behavior  of  HEMTs,  with  regards  to  signal 
distortion  and  intermodulation.®^ 

Situations  of  low  nao  coupled  with  the  possibility  of  strain-fields  (and 
concomitant  local  shifts  in  the  energy  band),  can  lead  to  :  (i)  Increased  localization  of 
charge  within  the  2DEG,  (ii)  Spatially  non-uniform  screening.,  and  (iii)  Scattering 
enhancement  of  the  2DEG.  Under  such  a  scenario,  electrons  would  be  bunched  within 
pockets  of  energy  minima.  Electron  flow  would  then  be  affected  by  the  Coulomb 
repulsion  from  charge  localized  at  dislocations,  leading  to  disruption  in  directed  motion 
and  lower  mobilities.  A  rigorous  treatment  of  such  transport  conditions  would  need  the 
inclusion  of  quantization  effects  within  the  regions  of  localized  strain,  and  use  of  a 
Molecular  Dynamics  (MD)  approach^®.  However,  neither  has  been  applied  here  in  our 
calculations,  for  the  following  reasons.  The  maximum  shift  *V  in  the  conduction  band  is 
given  by:  *V  ~  acb(2Br)’^[(l-20/(l-0]>  where  b  is  the  Burgers  vector,  ( the  Poisson’s 
ratio,  and  ac  the  band  deformation  potential.  Using  typical  values,  b  =  3.189  A  ,  (  =  0.3, 
and  ac  =  -8  eV  ,  leads  to  a  maximum  radial  distance  rmax  =  8.9  nm,  at  which  the  *V 
magnitude  equals  the  room  temperature  thermal  energy.  Hence,  it  is  expected  that  carrier 
at  distances  larger  than  rmax  would  not  be  localized  at  the  dislocations,  and  would  remain 
largely  unaffected.  Since,  this  rmax  value  is  much  smaller  than  the  typical  length  scales 
(e.g.  either  the  adjacent  dislocation  distances  of  ~  300  nm,  or  the  typical  HEMT 
dimensions),  such  quantization  effects  and  perturbations  in  local  density  were  ignored. 
Our  reason  for  not  using  the  MD  approach  was  to  avoid  computational  complexity. 
Besides,  MD  has  an  advantage  primarily  in  situations  that  involve  high  electron  and/or 
dislocation  concentrations.  However,  the  focus  here  is  on  low  electron  densities  for 
which  screening  and  effective  mobilities  are  the  smallest.  It  is  also  assumed  that 
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collective  dislocation  effects  would  not  be  an  issue  due  to  improvements  in  processing 
capability. 

Dislocation  scattering  was  treated  based  on  deformation  potential  interaction  of 
electrons  from  the  strain  fields  surrounding  edge  dislocations.  The  effect  of  a  strain  in  the 
quantum  well  is  to  shift  the  conduction  and  valence  band  edges.  Such  shifts  in  the 
conduction  band  edge  have  been  calculated, and  the  perturbing  potential  *V  for  electron 
scattering  can  be  shown  to  be  : 

|*V(r,2)|  =  {  ac  b  /(2Br)}  [  (l-20/(l-()]  sin(2) ,  (9) 

where  2  the  polar  angle  with  respect  to  the  Burgers  vector  in  the  2D  plane.  A  simple 
result  for  the  dislocation  scattering  rate  can  be  obtained  if  the  electronic  states  in  the  GaN 
channel  are  assumed  to  be  strictly  two-dimensional  in  nature.  For  this  case,  the 
wavefunctions  Pi(r)  for  the  i*''  subband  are  expressed  as:  Pj  (r)  =  exp(ik.r)  with  r  and 
k  being  the  position  and  wavevectors  parallel  to  the  interface,  and  S  the  surface  area.  The 
matrix  element  My  =  <Pi  |*V(r,2)|Pi  >  then  is: 

<Pi  |*V(r,2)|Pi  >  =  hV(q,N)|  =  {ac  b  /(Sq)}  [  (1-20/(1-01  sin(N)  ,  (10) 

where  N  is  the  angle  between  the  Burgers  vector  and  q  the  wavevector  magnitude 
between  the  initial  state  (k,)  and  final  in-plane  state  (Iq ).  Thus,  q  =  |  kj  -  k|  |.  Use  of  the 
Fermi’s  Golden-rule  for  this  two-dimensional  case,  with  an  additional  screening  factor  of 
k(q)  [=  1+  3/q]  for  the  matrix  element  (10),  then  leads  to  the  following  energy-dependent 
scattering  rate  Sdis(E)  for  parabolic  bands: 

Sdis(E)={Ndisac^  b^[(l-20/(l-0]^}[4BE  E]''  ol'  du  [{u+0.5S3/(2m*E)®-^}^(l-uY-^  ,  (H) 

where  Ndis  is  the  dislocation  density,  m*  the  electron  effective  mass,  and  3  the  inverse 
screening  length.  In  the  above,  an  average  over  the  angles  N  [defined  in  (10)]  has  been 
set  to  1/2  on  the  assumption  that  the  Burgers  vectors  associated  with  the  various 
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dislocations  are  randomly  oriented.  For  purposes  of  a  Monte  Carlo  implementation,  the 
final  in-plane  angle  after  a  dislocation  scattering  event  needs  to  be  determined.  This  is 
easily  done  based  on  the  associated  angular  probability  distribution  p(2).  In  this  case, 
from  (11) : 

p(2)  =  [(2  /  S)  {2m*E}"-^  sin(2/2)  +  3  .  (12) 

Hence,  for  the  j'*’  scattering  event,  the  requisite  scattering  angle  2j  can  be  obtained  by 
generating  a  random  number  q  ( 0  <  q  <  1),  and  using  the  above  p(2)  as: 

rj  =  [  oI  p(2)  d2  ]  /  [  oI  p(2)  d2  ]  .  (13) 

Due  to  the  strict  2D  nature  of  the  carrier  confinement  assumed  for  the  above, 
analysis,  the  scattering  rate  is  independent  of  the  sub-band  index.  A  more  accurate  result 
can  be  obtained  by  allowing  for  the  inclusion  of  the  envelope  function  along  the 
perpendicular  z-direction.  Using  the  following  separable  form^^:  Pi  (r)  =  exp(ik.r) 
Fi(z),  leads  to: 

|*V(q,N)|  =  {acb/S}[(l-20/(l-0]sin(N)ol"Fj*(z)  Fi(z)  dz  ol"  rJi(qr)/[r\zY-'  dr  ,  (14) 

where  J(x)  is  the  Bessel  function.  Integration  over  dr  yields: 

|*V(q,N)|  ={acb/S}[(l-20/(l-0]sin(N)ol"Fj*(z)  Fi(z)  x  [Io(x)  Ki(x)-Ii(x)  Ko(x)]  dz,  (15) 
where  x  =  qz/2,  the  Ii(x)  and  Ki(x)  are  Bessels  functions  of  imaginary  argument.  For  the 

92 

two  lowest  subbands,  the  following  expressions  for  the  envelope  wavefunctions  result  : 

Fi(z)  ~  [bo^/2]°-^  z  exp(-boz/2)  ,  (16a) 

and,  F2(z)  ~  [3biV{2(bo^  -bo  bj  [z  -  {(bo+  bi)/6}  z"  ]  exp(-biz/2)  ,  (16b) 
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with  fixed  parameters  bo  and  bi  that  depend  on  the  doping  and  channel  densities,  barrier 
thickness  etc..  As  such,  the  integral  in  (16)  can  only  be  evaluated  through  a  numerical 
procedure.  However,  useful  analytical  expressions  can  be  obtained  by  recognizing  that 
the  functions  Fj(z)  fall  off  rapidly  with  the  distance,  and  hence,  the  Bessel  functions  Ii(x) 
and  Ki(x)  can  be  approximated  for  small  x.  Using  simple  polynomial  expressions  for 
Ii(x)  and  Ki(x)^^,  and  employing  the  result  for  the  matrix  element  in  the  evaluation  of  the 
Fermi  golden-rule  scattering  rate,  leads  to  the  following  intra-subband  dislocation 
scattering  result : 

S”dis(E)~  Ndisac^  b^[(l-20/(l-0]^[4BSE]-^  ol^u  F[{u+.5S3/(2m*E)®-^}(l-uy^]-Mu  .  (17) 

Eqn.  (17)  is  almost  identical  to  (11),  except  for  the  additional  energy  dependent  factor 
F(E)  given  as:  F(E)  =  (1+  84.3  m*Eu^/{Zbo}^).  For  small  energies,  the  2D  and  quasi-2D 
scattering  rate  expressions  match.  The  increase  in  scattering  with  energy  is  associated 
with  the  larger  changes  in  wavevector  q.  The  angular  probability  function,  required  for 
selecting  the  final  state  after  dislocation  scattering  in  a  Monte  Carlo  implementation,  is 
modified  because  of  the  F(E)  factor.  Instead  of  (12),  the  angular  function  becomes; 

p”(2)  =  {1+  84.3  m*E  sin^(2/2)/  (Zbo)^}  [(2  /  Z){2m*E}®-^  sin(2/2)  +  3  f  .  (18) 

The  intra-subband  scattering  rate  S^^ciis(E)  can  also  be  obtained  by  using  the  F2(z) 
wavefunctions.  The  rate  for  the  intra-subband  transition  2->2  is  similar  in  form  to 
S^’dis(E),  except  for  a  different  F(E)  factor.  Denoting  this  factors  as  F^^(E),  one  obtains: 

F^^(E)  ~  A  (bo/bi  -  36.1468  m*E  u^/  {Zbi}^  )  ,  (19) 

where  A  =  [9b//{(bo^  -bo  bi  +bi^)^}].  The  angular  function  p^^(2)  is: 

p^^(2)={bo/bil  -  36.1468  m*E  sin^(2/2)/  (Zbi)^}  [(2  /  Z){2m*E}°-^  sin(2/2)  +3]'^.  (20) 
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The  inter-subband  rates  S^^dis(E)  and  S^^dis(E)  are  equal,  but  very  small  in  comparison  due 
to  minimal  overlap.  Hence,  these  have  been  ignored  here.  In  any  case,  the  dislocation 
scattering  rate  depends  on  the  bo  and  bj  parameters.  These  are  influenced  by  the  electric 
field  at  the  GaN-AlGaN  interface.  Hence,  it  becomes  clear  that  the  dislocation  scattering 
contributions  to  carrier  mobility  should  exhibit  a  dependence  on  both  the  Aluminum  mole 
fraction  and  barrier  thickness  as  these  affect  the  interface  electric  field. 

Computations  of  the  intra-  and  inter-subband  scattering  rates  require  evaluations 
of  the  bo  and  bi  parameters.  Approximate  solutions  can  be  obtained  on  the  basis  of  the 
energy  minimization  principle.^^  However,  a  slightly  different  approach  was  used  here. 
A  numerical  procedure  based  on  a  self-consistent  solution  of  the  Schr  'odinger,  Poisson, 
and  charge  balance  equations  was  employed  to  obtain  the  eigen-functions  explicitly. 
The  resulting  wavefunctions  were  then  used  to  extract  the  bo  and  bi  parameters.  Details 
of  this  method  for  the  GaN  electronic  system,  with  more  recent  inclusions  of  the  hole 
states,  have  been  reported  elsewhere.®"^'^^  Briefly,  both  spontaneous  and  piezoelectric 
polarizations  and  exchange-correlation  corrections  to  the  Coulomb  interaction  are 
included  in  the  calculation  of  the  electrostatic  potential.  Strain  is  included  in  the 
eigenstate  calculation  via  deformation  potential  theory.  The  large  band  gap  of  the  host 
materials  enables  the  conduction  band  to  be  decoupled  from  the  99v  and  two  37v 
valence  bands.  The  result  is  that  the  standard  8  x  8  k  •  p  Hamiltonian  inms  -p  basis  is 
transformed  into  two  equivalent  1x1  Hamiltonians  for  the  3jc  conduction  bands  and  a  6  x 
6  Hamiltonian  for  the  A,  B,  and  C  valence  bands.  The  electron  eigenstates  are  calculated 
by  solving  the  one-band  Schr'odinger  equation  across  a  c-axis  slice  of  the  HFET.  The 
hole  eigenstates  are  calculated  along  the  same  axis  using  a  6  x  6  k  •  p  method.  The 
complete  electron  distribution  was  calculated  by  summing  the  individual  contributions 
from  about  25  electron  subbands.  Summations  for  the  hole  distribution  were  performed 
over  70  hole  subbands  each  containing  an  admixture  of  Sqv  and  37v  states.  Also,  surface 
states  were  included  using  Fermi-Dirac  statistics,  since  surface  donors  strongly  influence 
the  2DEG  density  and  explains  the  experimental  data  quite  well.  The  resulting  model 
allows  for  the  calculation  of  electron  eigenfunctions,  energies  and  the  free  electron 
distributions  as  a  function  of  the  AlGaN  barrier  thickness  and  alloy  composition. 
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Fig.  4  shows  the  scattering  rate  as  a  function  of  energy  with  3  as  the  variable 
parameter.  For  calculations,  the  dislocation  density  was  2  x  10^  cm'^.  Since  the  scattering 


Fig.  4.  Energy  dependent  dislocation  scattering  at  300  K  with  the  inverse  screening 
length  as  the  parameter.  A  dislocation  density  of  2  x  10^  cm‘^  was  used. 

is  largest  at  small  energies,  one  can  expect  the  impact  to  be  strongest  on  carrier  mobility, 
rather  than  the  high-field  drift  velocity.  Also,  the  effects  would  be  greater  at  low  qs,  and 
hence,  for  situations  in  which  the  average  density  of  mobile  electrons  in  the  channel  was 
low.  High  defect  densities  (which  might  trap  free  electrons),  or  application  of  a  positive 
gate  bias  both  provide  such  conditions.  For  HEMTs  with  a  capping  layer  above  the 
AlGaN  barrier,  a  positive  bias  would  drive  holes  from  the  capping  layer-AlGaN  interface 
towards  the  2DEG.  For  situations  in  which  the  capping  layer  was  absent,  the  positive 
voltage  would  lower  the  surface  states  in  energy  space,  increase  their  occupancy  and 
thereby,  lower  nao  • 

4.  INTERMODULATION  DISTORTION  ANALYSIS 

Mathematical  evaluation  of  the  nonlinear  characteristics  of  field  effect  transistors 
(FETs)  has  often  been  based  on  the  Volterra-Wiener  formalism.^^’^^  However,  this 
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technique  is  only  applicable  to  small-signal  inputs,  and  it  neglects  memory-effects.  A 
more  appropriate  scheme  for  an  analytical  solution  would  be  a  polynomial  representation 
of  the  device  output  current  (or  current  density  Jo  )  as  a  function  of  the  input  gate  voltage 
Vg.^^  Thus,  Jo(Vg)  =  Ko  +  KjVg  +  KzVg^  +  KsVg^  +  ...+  with  the  coefficients  Ki  being 
some  constants.  For  accurate  and  realistic  results,  the  coefficient  can  either  be  extracted 
from  actual  experimental  data,  or  from  curve-fitting  to  Monte  Carlo  predictions  of  the 
current  response.  For  a  single  sinusoidal  input,  VG(t)  =  A  cos(Tt),  the  polynomial 
representation  yields  the  following  current  density  function  to  fifth  order: 

Jo(A,t)  =  (1/2)K2A^  +  (3/8)  K4A^  +  {KjA  +  (3/4)K3A^  +  (5/8)K5A^}  cos(Tt)  + 

+  {(1/2)K2A^  +  (1/2)K4A^}  cos(2Tt)  +  {(1/4)K3A^  +  (5/16)K5A^}  cos(3Tt)  +  (21) 

+  {(1/8)K4A^  }  cos(4Tt)  +  {(1/16)K5A^  }  cos(5Tt)  . 

For  a  two-tone  input  gate  signal  VG(t)  =  A  cos(Tit)  +  A  cos(T2t)  ,  the  current  density 
terms  at  the  fundamental  and  intermediate  frequencies  based  on  a  fifth  order  polynomial 
are: 


Jo(A,t)  =  {Ki  +  (9/4)  K3A^  +  (25/4)  KsA^}  [cos(Tit)  +  cos(T2t)]  + 

+  {(3/4)K3A^  +  (25/8)K5A^}  [cos{(2Ti-  T2  )  t}  +  cos{(2T2  -  Ti )  t}]  +  . . .  +  .  (22) 

The  choice  to  fifth  order  is  meant  to  represent  an  acceptable  tradeoff  between  increasing 
accuracy  with  the  inclusion  of  higher  orders  at  the  expense  of  enhanced  computational 
complexity.  Based  on  the  above,  expressions  for  the  input  power  (Pi),  the  output  power 
(Po),  in  dBm  units,  at  the  fundamental  frequency  can  be  obtained.  These  are: 

Pi  =  10  logio  [A^  10^/(2Ri)]  -30  +  10  logio  [A^  ]  ,  (23) 

and, 

Po  =  10  logic  [{KiA  +  (3/4)K3A^  +  (5/8)K5A^}^  10V(2Ro  )] 

-30+10  logio  [{KiA  +  (3/4)K3A^  +  (5/8)K5A^  f  ]  ,  (24) 
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where  Rf  and  Ro  are  the  respective  input  and  output  resistances.  Similarly,  the  variation 
of  the  third-order  intermodulation  power  Pip  ,  in  dBm,  with  the  gate  signal  amplitude  “A” 
is : 


Pip  -30  +  10  logio  [{(3/4)K3A^  +  (25/8)K5A^  f  ]  .  (25) 

Based  on  the  above  equations,  the  input-output  power  curves  can  be  obtained,  provided 
the  coefficients  Ki  are  given  or  can  be  extracted. 

The  above  approach,  though  reasonable,  excludes  memory  effects  within  the  system. 
It  is  implicit  that  the  output  changes  instantaneously  in  response  to  the  fluctuations  of  the 
input  signal.  A  more  realistic  representation  is  provided  through  Monte  Carlo 
simulations  that  can  yield  an  output  current  due  to  a  collective  motion  of  the  2DEG 
ensemble  in  response  to  gate  voltage  variations.  The  history  of  individual  particles  in 
automatically  included,  as  are  all  details  of  scattering  that  control  the  dynamical 
evolution  of  the  carriers.  The  non-Markovian  behavior  is  expected  to  become  important 
for  short  devices  and/or  high  operating  frequencies.  Here,  both  the  analytical  results  and 
the  Monte  Carlo  calculations  have  been  used  to  probe  the  IMD  behavior. 

5.  ANALYTICAL  MODEL  FOR  INTERMODULATION  DISTORTION 

A  more  appropriate  and  simpler  scheme  for  evaluation  would  be  through  an 
analytical  solution  would  be  a  polynomial  representation  of  the  device  output  current  (or 
current  density  Jo  )  as  a  function  of  the  input  gate  voltage  Thus,  Jo(Vd,Vg)  = 

ns(VD,  Vg)  q  v(Vd  Vg)  ,  where  ns  is  the  inversion  charge  density  and  v(Vd  Vg)  the 
effective  drift  velocity. 

For  all  practical  purposes,  the  density  of  inversion  electrons  in  the  GaN  channel  is 
well  above  10^®  m'^  level.  This  has  also  been  verified  through  the  calculations  presented 
later  in  this  report.  At  these  high  densities,  equation  (7d)  can  roughly  be  linearized  and 
expressed  in  the  form:  Ef(x)  =  A(x)  ns(x).  Combining  equations  (7)  and  using  the 
linearized  expression  for  Ef(x),  yields  the  following  inversion  charge  density  relation: 

n,(x,y)  =  ^s(x)  -  [,o,(-^)/d]  {Nb(x)  +  Ef(x)  -  )Ec(x)  -  Vg  +  V(y)  }  ,  (26a) 
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i.e.  ns(x,y)  =  Bi(x,Vg)  -  B2(x)  V(y)  -  [,o,W/d]  A(x)  n/x) ,  (26b) 

where,  Bi(x,Vg)  =  Bio(x)  +  [,o,(Jc)/d]  Vg  ;  B2(x)  =  [,o,(ji:)/d]  (26c) 

and  ,  Bio(x)  =  -  [,o,(-^yd]  {Nb(x)  -  )Ec(x)  }  .  (26d) 


Hence,  more  compactly : 

«s(x,y)  =  [Bio(x)  +  {,o,(x)/d}VG  -  B2(x)  V(y)]  /{I  +  [,o,(.^)/d]  A(x)  }  ,  (26e) 

which  can  be  used  in  the  evaluation  of  device  current. 

The  drift  velocity  term  v(Vd  Vg)  is  more  difficult  to  model  analytically.  Perhaps,  the 
best  approach  to  an  analytical  expression  for  the  drift  velocity  is  to  first  obtain  velocity- 
field  characteristics  based  on  a  rigorous  Monte  Carlo-type  simulation.  These  curves  can 
then  be  fit  to  polynomials  with  Vd  and  Vg  as  the  independent  parameters.  Unfortunately, 
the  use  of  such  polynomials  precludes  a  closed-form  analytical  expression  for  the  device 
current-voltage  characteristics.  For  example,  the  following  functional  form  has 
commonly  been  employed  for  the  velocity-field  (v-E)  curve: 

v(E)  =  [  pE  +  Vsa,  (E/Ec)"  ]  /  [1  +  (E/Ec)"]  ,  (27a) 

where  p  is  the  low-field  mobility,  Vsat  the  saturation  velocity,  Ec  the  electric  field 
corresponding  to  the  peak  velocity,  and  n  an  exponent  usually  equaling  about  4.  The 
above  form  includes  both  a  velocity  overshoot  behavior  and  high-field  saturation.  Using 
(26e)  and  (27a)  for  the  charge  density  and  velocity,  respectively,  then  yields  the 
following  expression  for  the  source-drain  current  Ids  • 

Ids  =  Zq[{Bi(x,VG)-B2(x)  V(y)}/{l+[,o,(x)/d]  A(x)}]  [pE  +Vsa,  (E/Ee)"  ]  /[1+(E/Ec)"] , 

(27b) 

where  Z  is  the  device  transverse  dimension,  and  q  the  electronic  charge.  Using  E  = 
dV(y)/dy,  equation  (27b)  can  be  set  up  as  an  integral  equation.  However,  an  analytical 
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solution  cannot  be  obtained,  and  so  a  closed-form  current-voltage  (Ids  -Vd)  relationship 
in  terms  of  Vo  does  not  result. 

This  difficulty  can  be  overcome  by  using  a  slightly  simpler  (and  more  approximate) 
expression  for  the  v-E  curve,  as  first  suggested  by  Trofimenkoff  A  closed  form 
relation  can  then  be  obtained.  The  Trofimenkoff  expression  is  ; 

v(E)=  pE/[l  +  (p/vsa,)E]  ,  (28) 

where  p  is  the  low-field  mobility,  and  Vsat  the  saturation  velocity.  Though  equation  (28) 
does  successfully  represent  a  low-field  mobility  and  a  high-field  saturating  behavior,  it 
does  not  mimic  the  velocity  overshoot  characteristic.  Hence,  strictly  the  Trofimenkoff 
expression  is  mostly  directly  applicable  to  semiconductors  such  as  Silicon.  In  the  present 
context  of  AIGaN-GaN  HEMTs,  the  above  velocity-field  model  can  perhaps  still  be  used 
as  a  simple  and  adequate  representation,  at  least  at  low  drain  voltages.  Furthermore, 
previous  Monte  Carlo  simulations  for  the  GaN  HEMTs  have  shown  an  absence  of 
velocity  overshoot'*’  .  Physically,  this  can  be  the  result  of  real  space  transfer  of  electrons 
from  the  GaN  channel  into  the  bulk  AlGaN  region.  The  higher  electronic  mass  in  AlGaN 
reduces  the  effective  electronic  drift  velocity,  and  effectively  suppresses  the  overshoot 
phenomena.  A  necessary  condition  for  this  is  that  the  T-L  valley  separation  be  greater 
than  the  GaN-AlGaN  conduction  band  offset  energy.  This  condition  is  easily  met  for  the 
GaN  HEMT  system  at  the  mole  fractions  typically  used,  and  so  the  velocity  overshoot  is 
strongly  curtailed.  Consequently,  the  use  of  the  simpler  Trofimenkoff  velocity-field 
expression  becomes  justifiable. 

Using  equation  (28),  the  current  equation  can  be  cast  in  the  following  form: 

Ids  =Zq  [{Bi(x,Vg)-B2(x)  V(y)}/{l+[,o,(A:)/d]A(x)}]  [{pdV(y)/dy}/{l-h(p/vsa,)  dV(y)/dy}] 

■  (29) 

Integrating  the  above  between  y  =  0  and  y  =  L,  corresponding  to  voltages  V(0)  =  0  and 
V(L)  =  Vd  ,  one  obtains  the  following  current- voltage  characteristics; 


26 


Ids  =  {II  Z  q  Vd  [Bio(x)  +  {,o,(x)/d}VG  -  0.5  B2(x)  Vd  ]  }/  [L  +  (iiMat)  Vd  ]  •  (30) 


In  reality,  the  velocity  v(E)  is  not  controlled  by  the  channel  electric  field  E  alone,  as 
suggested  by  eqn.  (28).  Instead,  it  is  also  influenced  by  the  gate  bias  Vo  for  a  variety  of 
reasons.  First,  Vq  controlled  2DEG  density  variations  affect  carrier  scattering  through 
changes  in  the  screening  of  polar  interactions  and  degeneracy  modification.'^’  Next,  with 
increasing  gate  bias,  interface  roughness  scattering  can  change  nonlinearly  with 
transverse  electric  field.  This  would  lower  the  velocity  with  increasing  Vq.  The  degree 
of  real  space  transfer  of  electrons  from  GaN  into  the  AlGaN  layer  is  also  dependent  on 
Vg  .  Such  a  gate  voltage  influence  on  the  channel  drift  velocity  can,  therefore,  be  better 
represented  by  modifying  the  v-E  expression  (28)  to  include  Vg  as  a  controlling 
parameter.  Consequently,  the  mobility  and  saturation  velocities  would  not  be  constants, 
but  rather  become  gate  voltage  dependent  fitting  parameters,  modifying  the  velocity-field 
relation  (28)  as: 


v(E,  Vg)  =  E  Ki{Vg}  /  [1  +  E  K2{Vg}  ]  .  (31) 

In  the  above,  Ki{Vg}  and  K2{Vg}  represent  Vg  dependent  mobility  and  an  inverse 
effective  field,  which  can,  in  principle,  be  extracted  from  a  curve  fitting  procedure  to  the 
results  of  Monte  Carlo  simulations  at  various  gate  voltages.  More  generally,  of  course, 
Ki,2{Vg}  /  Ki,2{Vg,T}  since  the  mobility  and  saturation  velocity  are  functions  of  the 
operating  temperature  T  as  well.  Use  of  (31)  instead  of  (28),  modifies  the  Ids  -  Vd 
characteristics  to  the  following  form: 

Ids  =  (  Ki{Vg}  Z  q  Vd  [Bi(x,Vg)  -  0.5  B2(x)  Vd  ]  )/  [L  +  K2  Vd  ]  .  (32) 

In  the  above,  the  Bi(x,Vg)  term  depends  linearly  on  gate  voltage.  The  non-linearity 
arises  from  the  (Ki{Vg}/[L  +  Vd  K2{Vg}])  factor  and  depends  on  the  device  length  L, 
the  drain  bias  Vd  ,  operating  temperature  and  Aluminum  mole  fraction.  Denoting  this 
factor,  for  convenience,  in  terms  of  a  polynomial  in  Vg  with  coefficients  Ci(VD)  as: 
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Ki{Vg}/[L  +  Vd  K2{ Vg}]  /  Co(Vd)  +  Ci(Vd)  Vg  +  C2(Vd)  Vg^  +  CsCVd)  Vg'  +  . . .  , 

(33) 

the  various  frequency  components  of  Ids  can  be  obtained  for  sinusoidal  gate  inputs.  For  a 
large  signal,  constant  amplitude,  gate  voltage  Vg  =  A  cos(‘^ct) ,  the  intrinsic  device  current 
at  the  fundamental  frequency  to  lowest  order  would  be: 

Ids  ~  AZqVD  [Ci  {Bio(x)  -  0.5  B2(x)Vd}  +  (Co+O.TSCzA^ )  B2(x)  ]  cosC^ct)  ,  (34a) 
i.e.  Ids  ('c)  /  I(x,Vd)  cos("ct)  .  (34b) 

Based  on  the  above,  the  input  power  (Pi)  and  output  power  (Po )  in  dBm  units  become:  Pi 
=  10  logio  [A^  10V(2Ri)]  and  Po  =  10  logic  [I^(x,Vd)  10^/(2Ro  )]  for  input  and  output 
resistance  Ri  and  Ro .  Thus,  depending  on  the  sign  of  the  coefficients  Co  ,  Ci  and  C2  , 
there  could  either  be  gain  enhancement  or  compression. 

Inter-modulation  distortion  characteristics  naturally  result  by  considering  a  two-tone 
gate  input  signal  VG(t)  =  A  cos(Tit)  +  A  cos(T2t).  Based  on  (23)  and  (24),  device  current 
at  the  intermediate  frequencies  |(2Ti,2  -  T2,i  |  to  lowest  order  is  given  as: 

Ids  =  0.75A’  ZqVo  [C3{Bio(x)-0.5B2(x)Vd}+C2  BjCx)  ][cos{(2T,-T2)t}+cos{(2T2-T,)t}] 

(35) 

The  corresponding  inter-modulation  power  Pjmd  is: 

Pimd  =  101og,o[(0.75A'ZqVD[C3{Bio(x)-0.5B2(x)VD}+C2B2(x)])"loV(2Ro)]  .  (36) 

The  third-order  intercept  point  (IP3)  can  be  obtained  from  (34)  and  (36)  since  it  is 
defined  to  be  the  power  level  at  which  the  linear  extrapolation  of  Po  equals  the  IMD 
power.  In  this  case,  P1P3  equals: 

PiP3  =  10  logic  [  {4/3}  (ZqVD)"  {Ti(x)'  /  T2(x)}  {10V(2Ro  )}]  ,  (37a) 
where,  Ti(x)  =  Ci{Bic(x)- 0.5  B2(x)  Vd}  +  CoB2(x)  ,  (37b) 
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and,  T2(x)  =  C3{Bio(x)-0.5B2(x)Vd}  +  C2B2(x)  .  (37c) 


6.  RESULTS  AND  DISCUSSION 

6.1  Transport  Calculations  For  Bulk  GaN 

Results  of  the  Monte  Carlo  simulation  for  bulk  wurtzite  GaN  at  300  K  are  given 
and  discussed  first.  Some  of  the  transport  parameters  for  GaN  material  have  not  explicitly 
been  measured.  Hence,  as  a  first  step,  the  10,000-electron  Monte  Carlo  code  was  used 
without  any  quantization  effects  to  simulate  the  transport  in  bulk  material.  The  primary 
objective  was  to  start  with  as  many  of  the  known  parameters  for  wurtzite  GaN  as 
possible,  and  fine  tune  other  values  to  achieve  a  reasonable  fit  with  available  reports  in 
the  literature.  Figure  5  shows  the  transient  drift  velocity  of  the  electronic  ensemble  in 
bulk  GaN  for  different  values  of  the  electric  field  parameter.  An  upperbound  of  400 
kV/cm  was  chosen  to  prevent  the  electronic  energies  from  becoming  excessively  large. 
This  precaution  was  necessary,  since  the  present  simulation  does  not  include  a  full 
bandstructure  calculation,  but  instead  relies  on  a  simple  nonparabolic  approach  to  the 
energy  bands.  A  velocity  overshoot  arising  from  intervalley  transfer  is  immediately  seen 
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Fig.  5.  Monte  Carlo  results  of  the  transient  electron  drift  velocity  at  300  K  in  wurtzite 
GaN  for  different  fields. 

from  the  curves,  similar  to  the  occurrence  for  GaAs.  The  threshold  field  for  this 
transferred  electron  effect  is  predicted  to  be  around  170  kV/cm.  For  electric  fields  higher 
than  this  threshold,  there  is  a  monotonic  decrease  in  the  steady-state  drift  velocity, 
although  a  much  higher  transient  overshoot  results.  At  400  kV/cm,  for  example,  a  peak 
velocity  of  about  6.5  x  10^  cm/s  is  predicted.  The  central  conclusion  is  that  GaN  material 
can  be  well  utilized  for  nonequilibrium  transport  and  would  yield  fast  response  for 
scaled-down  structures.  Steady-state  values  at  the  high  fields  are  also  larger  than  those 
obtainable  for  GaAs.  The  field-dependent  velocity-field  curves  for  the  material  at  two 
different  temperatures  are  shown  in  Fig.  6.  A  fairly  large  drift  velocity  of  1.4  x  10^  cm/s 


Fig.  6.  Monte  Carlo  results  of  the  electron  velocity-field  characteristics  at  300  and  600 
K  for  bulk  GaN. 

is  predicted  for  temperatures  as  high  as  600  K,  while  the  room  temperature  peak  is  nearly 
2.6  X  10^  cm/s,  which  translates  into  a  30-percent  increase  over  that  for  GaAs.  These 
results  also  indicate  that  high-field  drift  velocities  for  the  wurzite  phase  are  larger  than 
those  reported  for  zinc  blende  GaN.  Although  experimental  reports  of  the  low-field 
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mobility  demonstrate  higher  values  for  zinc  blende  due  to  a  smaller  ®L-valley  effective 
mass,  the  wurtzite  phase  would  have  an  advantage  for  high  field  applications.  This  is 
because  the  satellite  valley  for  the  wurtzite  phase  is  higher,  at  around  1.95  eV  versus  the 
1.45  eV  inter- valley  energy  separation  for  zinc  blende.  These  conclusions  are  similar  to 
those  reached  previously  by  Kolnik  et  al.^°°  Finally,  the  Monte  Carlo  results  for  the 
lowest  >valley  occupancy  as  a  function  of  time  are  shown  in  Fig.  7  for  various  electric 
fields.  Fractional  occupancies  of  90  percent  and  higher  are  predicted  for  applied  electric 


Fig.  7.  Monte  Carlo  predictions  of  the  time-dependent  >-valley  occupancies  at  different 
fields. 

fields  below  150  kV/cm.  There  is  a  relatively  sharp  transition  into  the  satellite  valleys  for 
fields  around  180  kV/cm.  Steady-state  occupancies  as  low  as  35  percent  are  predicted  at 
our  highest  field,  400  kV/cm.  These  results  are  also  fairly  close  to  those  obtained  by 
Kolnik  et  al.,  on  the  basis  of  a  more  sophisticated  Monte  Carlo  scheme  that  included  a 
full  bandstructure.  This  agreement  establishes  the  validity  of  the  transport  parameters 
used  in  the  present  model  and  underscores  its  utility  as  a  simple,  computational  efficient 
tool  for  analyzing  transport  for  field  below  400  kV/cm. 
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6.2  Monte  Carlo  Results  for  HEMTs 

Simulations  were  carried  out  at  300  K  for  the  GaN-AlGaN  HFET  structure.  The 
polarization  effect  is  known  to  be  dominant  in  this  material  system,  and  so  should 
influence  both  the  carrier  sheet  density  and  interface  electric  field.  Since  the  carrier 
density  controls  the  device  current  and  mitigates  scattering  through  enhanced  screening, 
it  is  critical  to  correctly  predict  the  mobile  charge  density  as  a  function  of  the  geometry 
and  operating  conditions.  The  interface  field  is  also  an  important  transport  parameter, 
since  it  shapes  the  electronic  wavefunctions,  thereby  influencing  the  interface  scattering 
process.  Consequently,  the  sheet  density  and  interface  fields  were  calculated  as  a 
function  of  the  gate  bias.  The  AlGaN  mole  fraction  and  thickness  were  the  variable 
parameters.  The  results  are  shown  in  Figs.  8  and  9  for  an  undoped  system.  From  both 


Fig.  8  Gate  voltage  dependence  of  the  electron  density  for  various  AlGaN  mole 
fractions.  The  thickness  was  30  run. 

figures,  it  is  evident  that  the  2D  electron  density  ns  is  large  and  almost  an  order  of 
magnitude  higher  in  comparison  to  values  reported  in  the  literature  for  the  GaAs 
heterosystems.  Our  values  match  the  recent  experimental  reports  for  the  AlGaN-GaN 
HFET.  Figure  8  for  30-nm  barriers,  shows  that  as  the  AlGaN  mole  fraction  is  increased. 
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the  sheet  density  (and  hence  also  the  interface  field)  is  enhanced.  This  is  due  to  the 
increased  values  in  both  the  spontaneous  and  piezoelectric  polarization  terms.  Increases 
in  the  band  offset  with  mole  fraction  also  contribute  to  a  higher  value. 


Fig.  9  Gate  voltage  dependence  of  electron  density  for  various  Alo.isGao.gsN  layer 
thicknesses. 

Also  as  expected,  the  density  increases  with  bias.  The  four  curves  of  Fig.  9  show  the 
effect  of  variations  in  the  AlGaN  barrier  thickness  on  the  bias-dependent  sheet  density. 
The  A1  mole  fraction  was  kept  fixed  at  0.15.  The  largest  barrier  thickness  provides  the 
most  mobile  carriers  in  the  channel.  Obviously  there  would  be  a  practical  limit  on  barrier 
thickness  increases  for  purposes  of  enhancing  the  channel  density.  With  a  very  thick 
AlGaN  layer,  the  device  would  no  longer  be  pseudomorphic  and  would  have  internal 
dislocations  and  defects. 

Next,  the  electron  channel  mobility  was  calculated  based  on  the  Monte  Carlo 
model.  Currently  there  are  no  measured  data  available  on  the  interface  roughness 
characteristics,  to  the  best  of  our  knowledge.  In  the  absence  of  such  data,  the  interface 
roughness  correlation  length  L,  and  the  rms  height  of  interface  roughness  fluctuations  ), 
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were  taken  as  adjustable  parameters.  It  was  implicitly  assumed  that  the  interface 
roughness  process  could  be  adequately  modeled  in  the  Born  approximation  with  the  use 
of  weak-perturbation  theory,  and  that  localization  effects  were  negligible.  Furthermore, 
the  presence  of  any  spatially  distributed  interface  polarization  charge  arising  from  the 
creation  of  AlGaN  and  GaN  islands  at  the  hetero-boundary  were  presumed  to  be 
adequately  included  by  the  L  and  )  parameters.  Since  the  mobility  has  to  be  computed  at 
zero  (or  very  low)  electric  fields,  this  parameter  was  evaluated  based  on  the  diffusion 
coefficient  D.  Use  of  a  drift  velocity  calculation  at  low  fields  for  determining  the 
mobility  can  lead  to  large  statistical  variations,  due  to  the  low  driving  field.  The 
diffusion  constant  was  computed  from  the  time  derivative  of  the  second  central  moment 
as : 


£)=  0.5d[<(xi(ij- <x,  (2‘)>)^>]/dt  =  dP(fj/dt  .  (38) 

Strictly,  the  above  applies  to  nondegenerate  situations.  For  a  more  precise  calculation  of 
D  with  degeneracy,  a  modified  Monte  Carlo  technique  as  suggested  by  Thobel  et  al.'°’ 
needs  to  be  used.  However,  this  aspect  was  ignored  here,  and  instead  equation  (38)  used 
for  calculations  of  D.  An  approach  similar  to  the  present  one  was  taken  by  Yamada  et 
al.  in  the  recent  past.  The  low-field  mobility  .•  is  related  to  D  via  the  Einstein  relation 
as  D  =  (n:/q)  [dn/dE/r]'^  where  q  is  the  electronic  charge,  n  the  carrier  density,  and  Ep  the 
Fermi  energy.  For  a  two-subband  electron  system  considered  here,  the  density  n  is 
related  to  the  Fermi  level  as  : 

n  =  [m  ksT  /(BS^)]  ln[(l  +  exp{(EF-Ei)/(kBT)})  .  (1  +  cxp{(EF-E2)/(kBT)})]  ,  (39) 

where  m*  is  the  effective  electron  mass,  while  Ei  and  E2  are  the  subband  energies.  With 
the  use  of  the  above  equations,  the  mobility  can  be  determined  from  a  knowledge  of  the 
diffusion  constant  D.  Figure  10  shows  the  time-dependent  variation  of  the  P(t)  term  for 
various  values  of  L  and  )  with  applied  gate  bias  values  of  0.0  and  3.0  V.  A  30-nm 
Alo.15Gao.85N  layer  was  used.  Values  of  1.5  nm  and  2  nm  for  L,  with  )  =  1  and  0.75  nm, 
have  been  shown.  These  parameters  were  adjusted  until  a  good  match  of  the  electron 
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mobility  with  available  experimental  data  was  obtained.  The  best  mobility  fit  was 
obtained  for  Z,  =  1.5  nm  and  )  =  0,65  nm.  The  mobility  results  at  300  K  as  a  function  of 


Fig.  10.  Time  dependence  of  the  second  central  moment  obtained  from  Monte  Carlo 
calculations. 
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Fig.  11.  Monte  Carlo  results  of  the  gate-bias-dependent  GaN  HFET  electron  mobility 
with  and  without  degeneracy. 

gate  bias  are  shown  in  Fig.  11,  along  with  experimental  data  taken  from  the  literature.'®^' 
The  measured  values  do  have  a  significant  spread. 

The  role  of  degeneracy  is  clearly  brought  out  in  Fig.  11.  Without  degeneracy,  the 
mobility  is  predicted  to  be  much  lower,  since  the  possibility  of  final  state  occupancy  is 
then  discounted.  The  predicted  values  agree  with  the  available  experimental  data,  but 
only  when  the  Pauli  exclusion  is  taken  into  account.  It  may  be  mentioned  that  the 
experiments  measured  Hall  mobility,  while  here  the  low-field  drift  mobility  has  been 
obtained.  The  two  differ  by  the  Hall  factor.  The  Hall  factor  can  be  quite  complicated, 
since  it  depends  on  the  details  of  the  prevalent  scattering  mechanisms,  their  relative 
strengths,  and  the  bandstructure  that  the  mobile  carriers.  Typically,  the  Hall  factor  has  a 
value  higher  than  unity.  In  some  cases,  however,  such  as  holes  in  strained  Si  or  Sii.;cGex 
alloys,  this  factor  can  be  below  unity. This  is  due  to  the  strong  nonparabolicity  and 
band  warping  effects.  For  electron  in  GaN,  however,  such  a  peculiar  situation  does  not 
arise.  Hence,  the  Hall  factor  is,  in  general,  larger  than  1.  Consequently,  the  drift  mobility 
is  expected  to  be  lower  than  the  Hall  mobility.  This  aspect  was  indirectly  taken  into 
account  here,  and  the  fitting  parameters  were  chosen  to  yield  drift  mobility  predictions 
that  were  somewhat  lower  than  the  experimental  values. 

Next,  the  electron  drift  velocity  was  calculated  for  the  interface  roughness 
parameters  obtained  from  the  fitting  procedure  above.  The  transient  velocities  at 
different  values  of  longitudinal  electric  field  along  the  channel,  with  the  applied  gate  bias 
Vg  as  a  parameter,  are  shown  in  Fig.  12  for  300  K.  As  expected,  the  highest  velocity 
occurs  at  the  lowest  (0  V)  gate  bias,  due  to  reduced  interface  roughness  scattering.  At 
electric  fields  of  20  and  50  kV/cm,  the  steady-state  velocities  for  Vg  =  0  are  predicted  to 
be  about  1.5  x  lO’  and  2.1  x  10^  cm/s,  respectively.  At  the  3  V  gate  bias,  the 
corresponding  values  for  the  20  and  50  kV/cm  fields  are  1.15  x  10  and  1.8  x  10  cm/s. 
The  steady-state  results  gleaned  from  the  transient  Monte  Carlo  simulations  are  shown  in 
Fig.  13  for  four  values  of  the  gate  bias.  In  addition,  steady-state  velocity-field  curves  for 
the  2DEG  at  600  K  are  also  shown.  The  velocity  decreases  with  the  gate  bias  for  every 
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longitudinal  field  value,  but  the  scaling  is  not  linear.  This  is  expected,  since  changes  in 
Vg  primarily  affect  interface  scattering,  while  the  other  processes  are  only  mildly  changed 
due  to  variations  in  the  wavefunction  overall  and  screening.  In  any  case,  it  becomes  clear 


Fig.  12.  Transient  2D  electron  drift  velocities  as  a  function  of  the  gate  bias  and 
longitudinal  electric  field  at  300  K. 
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Fig.  13.  Steady-state  velocity-field  characteristics  for  2D  electrons  for  various  gate 
voltages.  Temperatures  of  300  and  600  K  were  used. 

that  the  nonlinear  dependence  on  gate  voltage  is  likely  to  play  a  role  in  large  signal 
amplification  and  cause  mixing  of  multifrequency  signals.  The  values  for  600  K  are 
lower  than  at  300  K,  but  the  nonlinearity  with  Vg  appears  to  be  somewhat  mitigated.  This 
is  to  be  expected,  since  the  phonon  processes  become  more  dominant  at  the  higher 
temperatures  and  reduce  the  relative  role  of  interface  roughness  scattering  (IRS).  Since 
the  IRS  rate  is  strongly  energy  dependent  and  decreases  sharply  with  increasing  energy, 
the  gate  voltage  effects  should  be  most  prominent  under  conditions  of  low  carrier 
energies,  such  as  with  low  longitudinal  electric  fields.  This  is  borne  out  in  Fig.  13.  At 
the  lowest  electric  field  of  10^  Vm'\  the  drift  velocities  V300  and  Vm  corresponding  to 
temperatures  of  300  and  600K,  respectively,  are  predicted  to  be  1.044  x  10'^  and  0.4883  x 
10'^  ms'^  for  Vg  =  0.  However,  at  Vg  =  4.5  Volts,  the  velocities  V300  and  V^oo  are  computed 
to  be  0.59  X  10^  and  0.3018  x  10^  ms‘\  respectively.  The  ratio  {V600IV300)  |vg=o  then  turns 
out  to  be  less  than  (ymiysoo)  lvg=4.5-  This  in  inequality  (y6oolV3oo)  !vg=o  <  iV60olV3oo) 
|vg=4.5  is  the  direct  consequence  of  a  disproportionately  large  reduction  of  the  low-field 
300  K  velocity  in  going  from  the  Vg  =  0  operating  condition  to  Vg  =  4.5  Volts. 
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6.3  Monte  Carlo  Based  Distortion  Analysis 

The  polarization  effect  is  known  to  be  dominant  in  the  GaN-AlxGai-xN  material 
system,  and  influences  the  carrier  sheet  density  ris,  screening  and  degeneracy,  the 
interface  electric  field,  and  the  interface  roughness  scattering.  The  calculated  carrier 
density  as  a  function  of  the  gate  bias  for  an  undoped  system  with  30  nm  AlxGai-xN 
barriers,  has  been  shown  in  Fig.  8.  It  shows  a  negative  gate  threshold  for  each  curve, 
making  its  a  strongly  nonlinear,  monotonically  increasing  function  of  Vo  •  For  analytical 
calculations,  polynomial  curve  fits  to  the  aIs-Vg  relationships  can  easily  be  obtained  for 
each  mole  fraction.  Obviously  if  the  HEMT  were  operated  in  the  large-signal  mode  with 
amplitudes  less  than  the  gate  threshold,  then  the  nonlinearity  (and  hence,  IMD)  would  be 
less  severe.  Figure  8  also  indicates  that  increasing  the  mole-fraction  should  improve  the 
IMD  characteristics,  and  permit  higher  signal  voltages  and  input  power  due  to  an 
increased  gate  threshold.  Obviously  there  is  a  practical  limit  on  barrier  thickness 
increases,  since  for  very  large  mole  fractions,  the  device  would  no  longer  be 
pseudomorphic,  and  could  have  internal  dislocations  and  defects. 

The  second  source  of  nonlinearity  arises  from  a  Vo  dependence  of  the  2DEG  drift 
velocity.  This  aspect  was  probed  by  carrying  out  Monte  Carlo  simulations  for  a  fixed 
drain  field  for  various  values  of  dc  gate  voltages.  Results  showing  variations  in  the 
electron  drift  velocity  Vdr  with  gate  voltage  are  shown  in  Fig.  14  with  the  mole  fraction  as 
a  parameter.  The  temperature  was  set  to  300  K  and  a  longitudinal  channel  field  of  10 
kV/cm  was  used.  The  relatively  low  field  was  chosen  to  avoid  complications  arising 


39 


3.5 


Gate  Voltage  (Volts) 

Fig.  14.  Monte  Carlo  results  of  the  electron  drift  velocity  at  300  K  as  a  function  of  the 
HEMT  gate  voltage  with  the  AlxGai-xN  mole  fraction  as  a  parameter. 

from  real  space  transfer  and  trapping  within  the  AlGaN  layer.  Such  secondary  effects 
have  been  reported  recently  at  high  drain  currents.^^  A  general  monotonic  decrease  of  the 
carrier  velocity  with  gate  voltage  is  apparent,  and  arises  from  increases  in  interface 
roughness  scattering  (IRS).  Since  the  decreasing  Vdr  behavior  is  opposite  to  the  its  trend, 
these  two  factors  mitigate  the  overall  nonlinearity  in  the  current.  Based  on  the  above 
argument,  one  can  expect  the  nonlinear  response  trend  to  change  at  higher  temperatures. 
This  is  because  at  higher  temperatures,  phonon  scattering  would  increase,  reducing  the 
role  of  IRS.  Consequently,  the  decrease  in  Vdr  with  gate  voltage  would  become  less 
significant,  and  the  Hs-Yg  term  would  dominate.  The  Vg  dependence  of  the  current 
density  was  obtained  by  combining  the  results  of  Figs.  8  and  14  for  various  values  of  the 
mole  fraction.  The  results  are  shown  in  Fig.  15.  The  x=0.15  curve  is  dominated  by  the 
electron  density  term.  It  has  a  strongly  non-linear  current  density  with  values  starting 
near  zero  close  to  the  -3.5  Volt  threshold.  This  is  followed  by  a  significant  increase.  The 
x=0.5  curve  allows  for  a  much  larger  gate  voltage  swing,  and  is  dominated  by  the 
velocity  term  since  carrier  density  does  not  change  much  over  this  range  of  Vq-  A 
monotonically  decreasing  trend  results.  Finally,  the  x=0.3  curve  exhibits  an  initial 
increase  due  to  the  /ij  term,  followed  by  a  decrease  arising  from  the  drift  velocity  factor. 
In  general,  the  results  exhibit  the  following  features,  (a)  All  curves  exhibit  non-linearity. 
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(b)  The  characteristics  change  with  the  mole  fraction  and  the  gate  voltage  swing.  For 
example,  the  x=0.3  curve  reveals  that  a  large-signal  gate  amplitude  limited  to  ±  2  Volts 
would  produce  a  nearly  constant  current  response.  For  these  conditions  then,  output 
power  at  the  fundamental  frequency  can  be  expected  to  be  relatively  low.  (c)  The  results 
of  Fig.  15  suggest  the  possibility  of  judiciously  selecting  conditions  for  minimal  non¬ 
linear  distortions. 


Fig.  15.  Current  density  dependence  on  the  gate  voltage  at  300  K  for  various  AlxGai-xN 
mole  fractions. 

The  input-output  power  characteristics  are  shown  in  Fig.  16.  Polynomial  fits  to  the 
voltage  dependent  current  densities  of  Fig.  15  were  used  to  generate  the  power 
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Fig.  16.  Calculated  output-input  power  characteristics  for  the  fundamental  and  third- 
order  inter-modulation  frequencies.  Mole  fractions  of  0.15,  0.30  and  0.50  were  used  for 
the  barrier  layer. 

characteristics  for  x=0.15,  0.30  and  0.50.  Values  of  Pin  on  the  x-axis  were  calculated  as: 
Pin  =  10  logio(10^  *A^)  with  A  being  the  gate  voltage  amplitude.  Similarly,  values  of  Po 
were  obtained  as:  Po  =  10  logio(10^  *J^)  from  the  computed  drain  current  density  J.  The 
following  features  of  interest  can  be  discerned  from  Fig.  16.  (i)  The  third-order 

intermodulation  curves  exhibit  a  local  minima  for  all  three  mole  fractions.  Such  behavior 
has  been  reported  for  GaAs  transistors^^  and  so,  is  consistent  with  experimental 
observations.  At  low  input  power,  the  curves  all  follow  the  simplified  theory  that  predicts 
a  constant  3  dB  increment  in  distortion  for  every  decibel  of  input  power,  (ii)  For  optimal 
results,  it  would  clearly  be  best  to  operate  the  device  at  an  input  power  level  that 
corresponds  to  the  IMD  minima,  since  the  differential  between  the  fundamental  and  IMD 
outputs  is  then  the  greatest.  This,  in  theory,  would  allow  for  superior  device  operation  at 
relatively  input  powers  larger  than  predicted  by  a  conventional  assessment,  (iii)  The 
lowest  power  output  at  the  fundamental  frequency  is  for  the  x=0.3  case.  This  outcome  is 
expected  from  Fig.  15,  which  as  discussed  previously,  reveals  a  nearly  constant  current 
response  for  large  signal  amplitudes  limited  to  2  Volts.  By  the  same  reasoning,  the 
fundamental  output  for  x=0.15  is  the  best,  since  the  corresponding  curve  of  Fig.  15  has 
the  most  linear  characteristic  of  the  three,  (iv)  Despite  low  output  power  at  the 
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fundamental  for  x=0.3,  there  would  seem  to  be  two  inherent  advantages.  First,  the 
characteristics  are  predicted  to  have  gain  expansion  at  higher  voltages.  Such  gain 
expansion  has  been  reported  in  the  literature^®^.  Second,  the  operating  power 
corresponding  to  the  IMD  minima  is  the  largest  for  x=0.3. 

The  large-signal  response  at  an  elevated  temperature  of  600  Kelvin  was  analyzed  next 
for  a  relative  comparison  with  room  temperature  results.  Simulation  predictions  obtained 
for  the  current  density  as  a  function  of  Vg  are  shown  in  Fig.  17  for  x=0.15.  As  expected, 
the  300  K  curve  (the  same  as  shown  in  Fig.  15)  is  generally  higher  due  to  increased 


Fig.  17.  Comparison  of  the  calculated  Vg  dependent  current  density  at  300  and  600 
Kelvin  for  an  Aluminum  mole  fraction  of  0.15. 

carrier  drift  velocity  resulting  from  lower  phonon  scattering,  except  near  the  Vg 
threshold.  Carrier  densities  at  the  two  temperatures  are  nearly  the  same,  expect  for  a 
slightly  more  negative  gate  threshold  for  the  600  K  case.  Hence,  the  curves  are  generally 
controlled  by  the  carrier  drift  term,  but  show  a  somewhat  larger  current  at  600  K  close  to 
the  gate  threshold  value.  Furthermore,  since  phonon  scattering  is  strongly  temperature 
dependent,  while  interface  roughness  scattering  is  not,  the  relative  role  of  IRS  is  reduced 
at  600  K.  Hence,  for  the  600  K  case,  the  velocity  variation  (and  thus  the  current  density) 
tends  to  flatten  out,  showing  a  weak  Vg  dependence  at  the  larger  gate  voltages.  Overall, 
the  600  K  characteristic  is  smoother,  and  is  perhaps  a  closer  approximation  to  a  linear 
curve  as  compared  to  the  300  K  result.  Thus  crudely,  one  might  expect  lower  distortion 
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for  higher  temperature  operation.  The  input-output  power  characteristics  obtained 
through  polynomial  curve  fits  to  the  data  of  Fig.  17  are  shown  in  Fig.  18.  Looking  only 
at  the  fundamental  frequency,  the  output  power  at  BOOK  exceeds  that  at  600K.  This  in 
itself  seemingly  favors  room  temperature  operation.  This  result  is  the  direct  consequence 
of  a  greater  drift  velocity  and  hence,  device  current,  at  BOOK  due  to  lower  total  scattering. 
However,  the  IPB  curve  for  600  K  is  substantially  smaller  than  the  corresponding  BOOK 
characteristic  and  exhibits  a  much  larger  differential.  Thus,  for  better  device  performance 
in  terms  of  an  IPB  differential  alone,  the  600  Kelvin  operation  seems  to  be  the  favored 
alternative.  Furthermore,  the  fundamental  and  IPB  curves  at  BOO  K  are  seen  to  intersect 
over  the  input  range  considered,  while  those  for  600K  do  not.  These  results  are  thus 
indicative  of  potential  benefits  for  device  operation  at  elevated  temperatures. 


Fig.  18.  Comparison  of  the  calculated  output-input  power  characteristics  for  the 
fundamental  and  third-order  intermodulation  frequencies  at  BOOK  and  600  K.  A  mole 
fraction  of  0.15  was  used. 

Finally,  direct  full-Monte  Carlo  calculations  for  the  HEMT  were  performed  to  gauge 
the  response  to  two-tone,  large-signal  inputs.  Instead  of  using  polynomial  fits  to  the 
steady-state  velocity-field  data  as  generated  by  a  Monte  Carlo  calculation,  a  full  dynamic 
MC  simulation  of  the  current  response  to  a  two-tone,  time  varying  signal  was  performed. 
The  MC-based  approach  was  undertaken  for  the  following  reasons,  (i)  It  provides  a 
much  better  dynamical  treatment  of  the  carrier  transport,  (ii)  The  results  serve  as  a 
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validity  check  to  the  analytical  fits,  (iii)  All  memory  effects,  which  were  absent  from  the 
polynomial  expressions,  could  be  included.  Also,  in  effect,  this  naturally  allowed  a 
voltage  dependence  to  the  various  coefficients  K;  .  (iv)  The  inclusion  of  distributed 
internal  effects.  In  the  analytical  approach  it  is  assumed  that  the  current  changes 
instantaneously  in  response  to  the  applied  signal  uniformly.  The  MC,  however,  takes 
account  of  the  spatial  distributions  and  temporal  delays.  The  results  for  the  fundamental 
and  IP3  are  shown  in  Fig.  19  at  300  K  for  x=0.15.  Running  the  MC  simulator  for  a 
discrete  set  of  Vg  large-signal  amplitudes  generated  the  required  data  points.  The  time- 
dependent  current  density  response  was  recorded  for  each  input  two-tone  signal.  Due  to 


Fig.  19.  Direct  Monte  Carlo  results  of  the  output-input  power  characteristics  at  300K  for 
a  mole  fraction  of  0.15. 

the  computationally  intensive  nature  of  the  simulations,  only  a  few  discrete  Vg 
amplitudes  were  chosen.  The  use  of  relatively  few  Vg  points,  coupled  with  the  statistical 
nature  of  the  MC  simulation,  yielded  curves  that  were  less  smooth  compared  to  those  of 
Fig.  16  that  resulted  from  the  analytical  fits.  Similarities  between  the  curves  of  Fig.  16 
and  the  300K  analytical  results  of  Fig.  18  are  apparent.  In  both  cases,  the  IMD  minima 
occurs  at  a  value  of  about  35,  though  the  analytical  curve  of  Fig.  16  is  much  smoother. 
Minor  differences  include  a  more  pronounced  saturation  (especially  at  the  fundamental) 
and  a  higher  cross-over  point  from  the  Monte  Carlo  calculations.  The  disparity  is 
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indicative  of  the  internal  memory  and  distributed  effects  that  only  the  MC  calculations 
can  fully  account. 

6.4  Analytical  Assessment  for  IP3  and  Distortion 

Equation  (37)  reveals  that  the  requirements  for  reduced  inter-modulation  distortion, 
increased  gain,  and  maximization  of  the  IPS  point  are  increases  in  Z,  Co ,  Ci  and  Ti(x), 
with  corresponding  reductions  in  T2(x).  The  latter  can  be  attained  either  by  designing  C2 
0  and  C3  ->  0,  or  more  realistically,  by  ensuring  that  C2  and  C3  have  opposite  sign. 
From  a  physical  device  standpoint,  possibilities  for  achieving  a  large  IPS  dynamic  range 
then  amount  to  the  following:  (i)  Reducing  the  channel  length  L  so  that  [Ki/{L+VdK2}] 
Ki/(VdK2)  VsatA^D  .  This  would  almost  be  a  fixed  value  (for  a  given  temperature) 
with  only  a  weak  dependence  on  Vd  .  (ii)  Attempting  to  increase  the  barrier  height  for 
reductions  in  real  space  transfer  (RST),  thereby  mitigating  the  role  of  gate  voltage  on  the 
transport  characteristics.  This  can  be  achieved  either  by  increasing  the  mole  fraction  or 
inserting  a  very  thin  AIN  interfacial  layer  as  recently  reported  by  Shen  et  al.^^°  The 
suppression  of  RST  would  effectively  make  Vsat  nearly  independent  of  Vg  ,  as  there  would 
only  be  one  conduction  channel  with  the  GaN  region,  (iii)  High  temperature  operation 
since  the  role  of  phonon  scattering  then  increases,  while  the  relative  contribution  of 
interface  roughness  scattering  (and  hence,  Vg)  is  reduced,  (iv)  Choosing  a  mole  fraction 
X  so  that  the  channel  carrier  density  does  not  reduce  to  zero  even  at  the  least  negative  Vg 
point.  On  the  basis  of  eqn.  (26),  this  implies:  Bio(x)  >  [{,o,(-^yd}|VG  |  +  0.5  B2(x)  Vd]. 
Mathematically,  this  translates  into  a  lower-bound  for  x  depending  on  the  maximum 
operating  bias  Vd  and  gate  voltage  swing  “A”  assigned  for  a  given  application.  Using  the 
mathematical  expressions  given  previously,  yields  the  following  inequality  for  improving 
the  dynamic  range: 

0.0483x{  [(5x  +  103)/(-32x  +  405)]  (0.73x+0.73  )  +  (0.1  Lc  +  0.49)  }+  0.052  x  > 
[8.85  X  10’^^  (9.5-0.5x)/d]  {0.84  +  O.lOSx  -  0.7x^  +  |Vg  |  +  0.5  Vds  }  ,  (40) 

Using  typical  upperbounds  of  10  Volts  for  Vds  and  2.5  for  Vg  ,  the  mole  fraction  needs  to 
be  greater  than  0.22  for  25  nm  AlGaN  thickness  and  above  0.49  for  a  10  nm  AlGaN 
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layer.  This  simple  calculation  helps  bring  out  an  inherent  advantage  of  the  GaN  HEMT 
for  microwave  applications.  As  is  evident,  the  existence  of  a  polarization  induced 
interface  charge  density  ^s(x),  enhances  Bio(x),  thereby  making  it  possible  to  avoid  a 
sharp  fall-off  in  the  channel  charge  for  more  negative  gate  voltages  even  under  large 
drain  bias  conditions.  There  is  thus,  a  significant  advantage  over  traditional  GaAs-based 
MESFETs  or  pseudomorphic  HEMTs. 

The  coefficients  Ci  were  extracted  from  Monte  Carlo  simulations  of  the  electron  drift 
velocity  characteristics.  For  purposes  of  evaluating  the  effects  of  real  space  transfer 
(RST)  on  devices  behavior,  Monte  Carlo  simulations  for  again  carried  out  at  300  K.  The 
corresponding  steady-state  velocity-field  curves  so  obtained  are  shown  in  Fig.  20.  For 
comparison,  results  without  RST  are  also  included.  The  RST  simulations  were  carried  out 
by  allowing  for  electron  transfers  into  and  out  of  the  AlGaN  barrier  layer.  The  AlGaN 
layer  was  treated  as  a  bulk  semiconductor,  and  quantization  effects  ignored  for  electronic 
transport  in  this  region.  The  velocity-field  (v-E)  curves  of  Fig.  20,  have  minimal 


Fig.  20.  Monte  Carlo  calculations  of  velocity-field  curves  with  inclusion  of  real  space 
transfer  at  300  K. 
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difference  at  Vg  =  0.  With  increasing  gate  voltage,  the  v-E  curves  are  predicted  to 
exhibit  gradual  reductions  in  slope  for  the  RST  situation.  This  is  consistent  with  the 
notion  of  redueed  “overall  mobility”  associated  with  parallel  conduction  in  both  bulk 
AlGaN  and  the  GaN  channel.  For  high  Vg,  however,  the  RST  velocities  are  somewhat 
larger  at  the  lowest  electric  fields.  This  is  indicative  of  lower  overall  scattering,  and  is 
associated  with  the  two  following  effects.  First,  interface  roughness  scattering  is  absent 
for  the  electronic  population  undergoing  bulk  AlGaN  transport.  Second,  at  low  electric 
fields  the  carrier  energies  are  lower  leading  to  the  reduced  probability  of  phonon 
emission  processes.  This  is  especially  relevant  for  this  material  system  given  the 
relatively  high  optical  phonon  energy. 

The  Monte  Carlo  velocity-field  curves  for  a  barrier  layer  mole  fraction  of  0.15  were 
used  to  determine  the  fitting  parameters  Ki{Vg}  and  K2{Vg}  of  eqn.  (31).  Results  of  the 
curve  fitting  at  300  K  and  600  K  with  and  without  the  inclusion  of  real  space  transfer 
(RST)  are  given  in  Tables  II(a)-II(c).  Based  on  the  data  presented  in  Tables  II(a)-(c),  it 
becomes  evident  that  the  variation  in  the  parameterized  ratio  R  (=Ki/K2)  is  the  least  for 
the  600  K  case,  and  the  most  for  at  300  K  with  inclusion  of  real  space  transfer.  For 
example  from  Table  11(c),  R{600K,  No  RST}has  a  very  modest  variation,  only  changing 
from  3.6  x  10^  at  Vg  =  0  to  3.67  x  10^  for  Vg  =  4.5  V.  The  deviation  in  R{300K,  RST} 
for  the  corresponding  gate  voltages  is  from  2.766  x  10^  to  1.946  x  10^  ,  while  R{300K, 
No  RST}ranges  from  2.6478  x  10^  to  3.475  x  10^  .  The  nearly  constant  value  of  R  at 
600K  obtained  here  bodes  well  for  reducing  gate-voltage  related  non-linear  distortions 
and  inter-modulation  effects.  So  at  least  for  small  channel  devices,  microwave 
performance  can  be  expected  to  improve  at  higher  operating  temperatures.  This  is  in 
keeping  with  the  previously  stated  expectation  of  a  relative  decrease  in  the  role  of 
interface-roughness  scattering.  The  results  at  300  K,  reveal  two  shortcomings  associated 
with  the  real  space  transfer  process.  First,  the  variation  in  R  with  Vg  is  higher  in  the 
presence  of  RST.  This  is  seen  more  clearly  from  the  following  parameterized  fit  for  R  at 
300  K  from  the  data: 

R{T=300,  No  RST}  =  2.64  x  10^  [1  -  0.142  Vg  +  0.0867  Vg^  -  0.00875  Vg^  ]  ,  (41a) 
and  R{T=300,  RST}  =  2.77  x  10^  [1  -  0.3646  Vg  +  0.1325  Vg^- 0.01469  Vg^]  •  (41b) 
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Table  11(a) 


Temperature  =  300  K ,  No  real  space  transfer ;  x  =  0.15 


Vg  (volts) 

0.0 

2.0 

3.0 

4.5 

Ki{Vg} 

1.88  X  10'^ 

1.23  X  10'^ 

0.944  X  10'^ 

0.695  X  10'^ 

K2{Vg} 

7.1  X  lO'"^ 

4.7  X  10'^ 

3.2  X  10’’ 

2.0  X  lO'"^ 

Table  11(b) 

Temperature  =  300  K ,  With  real  space  transfer ;  x  =  0.15 


Vg  (volts) 

0.0 

2.0 

3.0 

4.5 

Ki{Vg} 

1.66  X  10"^ 

1.77  X  10'* 

1.44  X  10'^ 

1.45  X  10'^ 

K2{Vg} 

6.0  X  lO'"^ 

9.4  X  lO'"^ 

7.4  X  lO'^ 

7.45  X  10'^ 

Table  11(c) 

Temperature  =  600  K ,  No  real  space  transfer ;  x  =  0.15 


Vg  (volts) 

0.0 

2.0 

3.0 

4.5 

Ki{Vg} 

0.648  X  10'^ 

0.484  X  lO'^ 

0.412  X  10'^ 

0.367  X  10'^ 

K2{Vg} 

1.8  X  lO'^ 

1.3  X  10’ 

0.9  X  lO  ’ 

1.0  X  10'’ 
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The  higher  order  non-linear  coefficients  are  obviously  larger  with  RST.  Second,  in  the 
presence  of  RST,  the  value  of  R  decreases  with  increasing  Vg  .  Such  a  trend  can  be 
expected  to  adversely  affect  the  gain  and  current  throughput  of  the  device  for  microwave 
signals  with  a  large  voltage  swing.  Thus,  it  seems  apparent  that  for  enhancing  the 
microwave  performance  of  GaN  HEMTs,  RST  should  be  minimized,  and  hence,  concepts 
such  as  the  inclusion  of  thin  AIN  barrier  layers  seem  very  beneficial. 

6.5  Transport  Calculations  For  GaN  Gunn  Structures 

Having  validated  the  Monte  Carlo  model,  simulations  were  then  carried  out  at  300 
K  for  the  GaN  based  Gunn  oscillatorbased  on  a  n^-n'-n-n"^  notch  structure  for  lengths  of 
the  n-region  around  1.0  :m.  The  much  smaller  length  (as  compared  to  the  3.0  :m  device 
simulations  previously  reported"^^  )  serves  two  purposes:  (a)  Enhanced  operating 
frequency  due  to  lower  transit  times,  and  (b)  ensures  a  higher  sinusoidal  component  at  the 
fundamental  mode.  Longer  devices  that  operate  in  dipole  domain  mode  tend  to  have  a 
flatter  response,  with  a  higher  non-sinusoidal  component  since  variations  in  terminal 
current  only  take  place  at  the  start  of  every  domain  formation  cycle.  The  n-doping  range 
of  5  X  10^^  m'^-  2  X  10^^  m'^  chosen  here,  leads  to  a  “nL”  product  above  the  5  x  10^^  m'^ 
threshold"^’  for  GaN  Gunn  oscillations.  For  a  more  realistic  simulation  of  the  microwave 
oscillator,  the  GaN  device  was  embedded  in  a  parallel  resonant  circuit  as  shown  in  Fig. 
21a.  This  allows  for  the  inclusion  of  circuit  effects,  permits  performance  evaluation  at 
the  fundamental  frequency  and  other  harmonics,  and  the  assessment  of  circuit  parameter 
variations.  Time-varying  voltage  boundary  conditions  were  applied  to  the  GaN  structure, 
and  the  device  currents  computed  based  on  the  Monte  Carlo  particle  flows.  The  total 
current  obtained  by  solving  the  related  Kirchhoff  equations.  The  displacement  current 
can  be  accounted  by  using  the  cold  capacitance  procedure^’ ^  In  our  simulations,  the 
device  structure  was  taken  to  have  a  cross-sectional  area  of  10’^  m'^  .  The  doping  profile 
for  the  conventional  notch-type  device  is  shown  in  Fig.  21b.  The  two  heavily  doped  end 
regions  with  densities  of  2  x  lO^'^  m‘^  serve  as  the  n^  contacts.  The  doping  level  “n”  of 
the  transit  region  was  taken  to  vary  within  the  5  x  10  m'  -  2  x  10  m'  range.  Thus, 
the  doping  level  of  the  drift  region  was  chosen  to  be  much  larger  than  typically  used  for 
GaAs  or  InP  Gunn  devices^^^’”^.  Since  the  mobility  of  GaN  is  lower  in  comparison  to 
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both  GaAs  and  InP,  the  time  required  to  build  up  a  stable  domain  is  higher  for  the  same 
background  doping.  The  higher  doping  chosen  here,  helps  offset  the  lower  mobility 
value.  The  notch  at  the  cathode  end  is  required  to  create  the  local  nonuniformity  that 
helps  launch  the  accumulation  or  dipole-domain  mode  instabilities.  The  equilibrium 
density  of  the  notch  was  set  at  20  percent  of  the  transit  region.  The  dead-zone  lies  in  this 
region.  Fig.  21c  shows  the  schematic  of  a  repetitive  Guim  diode  structure  with  three 
transit  regions.  Such  stmctures  with  multiple  zones  were  also  used  in  the  present 
simulations.  For  the  electrical  circuit,  a  20  Ohm  resistor,  an  external  capacitance  of  0.05 
pF  and  a  0.12  nH  inductor  were  assumed. 


Fig.  21.  Schematic  pertaining  to  the  device  simulations,  (a)  Lumped  electrical  circuit, 
(b)  Doping  profile  of  the  conventional  GaN  notch  diode,  and  (c)  Schematic  of  a  three 
transit  region  multiple  GaN  Gunn  diode. 
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The  circuit  of  Fig.  21  was  incorporated  and  the  initial  electron  distribution  taken  to 
mimic  the  doping  density  distribution.  The  Gunn  structure  was  taken  to  connect  with  the 
tank  circuit  at  t=0^  .  Monte  Carlo  results  at  300K  for  the  single  transit  Gunn  diode 
showing  the  time  dependent  behavior  of  the  total  current  and  the  device  voltage  are 
shown  in  Fig.  22.  The  device  length  was  taken  to  be  1.0  :m,  and  applied  bias  was  50 
Volts.  The  initial  transient  is  simply  associated  with  the  commencement  of  the  applied  dc 
bias  at  time  t  =  0.  A  steady  oscillatory  state  is  reached  in  about  60  picoseconds,  with  the 
operating  frequency  predicted  to  be  135  GHz.  Results  for  the  simulated  frequency 
dependencies  of  the  generated  power  and  device  efficiency  were  then  obtained  by 
running  such  time  dependent  simulations  for  various  device  voltages,  device  lengths  and 
operating  temperatures.  The  results  are  given  in  Figs.  23  and  24,  for  two  operating 
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Fig.  22.  Simulation  results  showing  the  temporal  development  of  the  device  voltage  and 
total  circuit  current  for  a  1.0  micron  device  at  300K.  A  50  Volt  bias  was  used. 

temperatures  of  300  K  and  450  K.  At  very  low  biasing  levels,  the  electric  field  within  the 
GaN  device  is  close  to  the  threshold  for  NDR,  and  hence,  the  power  output  is  low.  With 
increasing  bias  the  output  power  increases,  as  does  the  operating  frequency  due  to  the 
larger  average  electron  velocity.  A  central  maxima  thus  results,  and  the  general  shape  is 
in  keeping  with  previous  reports.  At  very  high  biasing  the  transferred  electron  effect  is 
lost  with  nearly  all  carriers  residing  primarily  in  the  satellite  valleys.  The  curve  at  the 
higher  temperature  of  450  K  shifts  to  lower  frequencies  due  to  the  reduced  carrier 
velocity  and  the  power  output  is  predicted  to  drop.  Fig.  23  shows  peak  powers  to  be 
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Fig.  23.  Monte  Carlo  results  of  the  power  output  versus  frequency  at  300K  and  450  K 
for  single  domain  GaN  Gunn  diodes. 

about  0.665  W  and  0.525  W  at  the  300  K  and  450K  operating  temperatures.  Since  device 
self-heating  effects  have  not  been  included  here,  these  values  represent  an 
overestimation.  In  any  case,  there  is  considerable  improvement  over  the  reported  results 
of  about  0.150  W  for  InP  This  is  not  surprising,  since  the  threshold  field  for  GaN  is 
roughly  ten  times  larger  than  in  InP.  The  corresponding  efficiencies  are  shown  in  Fig. 
24  at  the  300K  and  450  K  temperatures.  Peak  efficiencies  of  1.9  %  and  1.47  %  are 
predicted  for  the  single  transit  region  GaN  Gunn  diode.  The  room  temperature  value  is 
higher  than  that  reported  for  GaAs  devices^^^. 

Finally,  Monte  Carlo  simulations  were  performed  for  GaN  diodes  with  multiple 
transit  segments.  The  resulting  power  output  as  a  function  of  the  frequency  is  given  in 
Fig.  25  at  300  K  for  one,  two,  and  four  transit  regions.  A  simplistic  expectation  is  for  the 
output  power  to  scale  as  N^  (with  N  the  number  of  segments)  for  fixed  values  of  the 
operating  frequency,  temperature,  and  device  impedance.  This  follows  from  a  required 
scaling  in  both  the  voltage  (to  maintain  the  device  electric  fields)  and  device  current  (due 
to  an  increase  in  area  for  maintaining  the  impedance).  However,  in  practice,  the  benefits 
of  a  multiple  transit  structure  might  not  be  quite  as  significant  for  a  variety  of  reasons.  In 
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Fig.  24.  Simulation  results  of  the  device  efficiency  versus  frequency  at  300K  and  450  K. 

a  conventional  single-transit  region  device,  the  injected  carriers  are  all  nearly 
thermalized,  and  have  comparable  energies.  Hence,  the  device  noise  can  be  expected  to 
be  lower  due  to  a  relatively  “tightly-bunched”  carrier  distribution,  especially  for  short 
device  lengths.  With  multiple  transit  zones,  however,  the  carrier  distributions  at  the  start 


Fig.  25.  Results  of  the  power  output  versus  frequency  at  300K  for  multiple  transit 
segment  GaN  Gunn  diodes.  Plots  for  one,  two  and  four  segments  are  shown. 

of  each  successive  transit  segment  might  not  be  completely  periodic,  nor  closely  bunched 
in  distribution,  and  could  be  non-thermal.  Consequently,  one  can  expect  a  larger  variance 
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in  the  intervalley  transfer  distance,  leading  to  higher  noise.  The  actual  condition  at  the 
start  of  successive  domains  depends  on  the  thickness  and  doping  of  layer  between  the 
transit  regions.  Thicker  and  more  strongly  doped  layers  would  facilitate  better  carrier 
thermalization,  but  would  increase  the  overall  transit  time  and  adversely  impact  the 
operating  frequency.  The  absence  of  a  fully  thermal  distribution  at  the  beginning  of  each 
successive  transit  region,  also  implies  that  the  average  carrier  energies  within  the  device 
would  tend  to  be  larger.  This  could  lead  to  higher  internal  power  losses  because  of  strong 
polar-optical  phonon  scattering,  especially  within  the  satellite  valleys  due  to  their  higher 
density  of  states.  In  addition,  with  the  hot  electrons  spending  a  longer  fraction  of  their 
time  within  the  satellite  valleys,  the  overall  speed  would  be  slower,  and  the  swing  in  the 
current  oscillation  reduced.  Last,  but  not  least,  the  higher  power  losses  can  be  expected 
to  lead  to  higher  internal  temperatures  and  possible  thermal  gradients  due  to  local  hot 
spots  near  the  successive  “virtual  anodes”.  Both  would  adversely  affect  the  transport  and 
reduce  operating  frequency.  Collectively  then,  the  useful  power  output  might  not  scale 
quite  as  well  as  . 

Most  of  the  above  issues  could  be  addressed  by  a  three-dimensional  Monte  Carlo 
simulation  procedure.  However,  since  we  have  not  performed  self-consistent  electro¬ 
thermal  calculations,  the  heating  aspects  and  their  role  in  mitigating  the  Gunn  diode 
response  could  not  be  addressed.  The  thermal  analysis  as  well  as  the  noise  characteristics 
will  be  presented  elsewhere.  The  Monte  Carlo  results  of  the  power  output  versus 
frequency  are  given  in  Fig.  25  at  300  K  for  one,  two,  and  four  transit  regions.  The  device 
area,  however,  was  not  increased  for  ease  of  comparison.  The  results  show  that  the 
output  power  does  not  quite  scale  as  N,  but  instead  is  a  sub-linear  function.  Hence,  by 
increasing  the  area,  one  would  not  expect  the  output  power  enhancement  to  be  quite  as 
large  as  N^.  In  any  case,  the  power  output  is  certainly  better  than  that  of  the  conventional, 
single-transit  structure,  and  so  the  multiple  GaN  Gunn  diode  structure  does  merit  serious 
experimental  study. 

6.6  Transport  Effects  in  GaN  HEMTs  With  Dislocation  Scattering 

In  order  to  probe  dislocation  scattering  effects,  a  simple  HEMT  structure  consisting 
of  an  undoped  20  nm  Alo.1Gao.9N  barrier  on  GaN  was  assumed.  Surface  donor  states  of 
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density  5  x  cm'^  with  an  activation  energy  of  1.4  eV  were  assumed.  Fig.  26  shows 
the  resulting  zero-bias  energy  band  diagram  for  this  structure.  A  net  2DEG  of  1.6566  x 
lO’^  cm'^  resulted,  with  a  spatial  distribution  given  in  Fig.  27.  Due  to  the  finite  barrier 
height,  some  non-zero  electron  concentrations  in  the  Alo.1Gao.9N  layer  are  obvious.  A 
relatively  low  mole  fraction  was  used  to  reduce  the  lattice  mismatch  and  hence,  the  strain 
within  the  layers.  The  associated  electronic  wavefunctions  for  the  above  HEMT  structure 
for  the  first  three  subbands  is  shown  in  Fig.  28.  The  eigen-energies  of  these  three 
subbands  were  found  to  be  7.9068  meV,  57.34  meV,  and  81.238  meV,  respectively.  The 
parameters  bo  and  bi  given  in  eqn.  (9)  for  the  envelope  wavefunctions  were  then  obtained 
from  the  plots  of  Fig.  28,  through  a  simple  least-square  fitting  procedure.  Thus,  bo  ~ 
1.1428  X  10^  m'^  and  bi  ~  7.634  x  10®  m"^ . 

The  electron  mobility  in  the  HEMT  channel  for  the  above  device  consisting  of  an 
undoped  20  imi  Alo.1Gao.9N  barrier  on  GaN,  was  next  obtained  through  Monte  Carlo 


Fig.  26.  Energy-band  diagram  of  an  undoped,  unbiased  GaN-AloiGao.9N  HEMT  structure. 
The  electron  density  distribution  within  the  HEMT  is  also  structure  are  shown. 
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simulations.  The  results  as  a  function  of  lattice  temperature  are  shown  in  Fig.  29  with  the 
dislocation  density  Ndis  as  a  parameter.  As  might  be  expected,  mobilities  in  the  absence 
of  dislocations  are  the  highest.  Their  room  temperature  value  is  predicted  to  be  around 
1.711  X  10^  cm^A^s,  which  is  fairly  close  agreement  with  the  experimental  data  reported 
by  the  Santa  Barbara  group^®’  At  77  K,  the  highest  channel  mobility  was  obtained  at 
about  8.440  x  10^  cm^A^s.  This  is  somewhat  lower  than  the  experimental  report,  and  is 
probably  due  to  our  over-estimation  of  interface  roughness  scattering  (IRS).  At  the 
highest  dislocation  density  of  10^°  cm“^,  mobility  reduction  of  16.8  percent  and  8.6  per- 


Fig.  29.  Electron  mobility  versus  temperature  with  dislocation  density  as  a  parameter. 

cent  are  predicted  for  77  K  and  300  K,  respectively.  Given  the  sharp  drop  in  phonon 
scattering  at  the  lower  temperature,  the  impact  of  dislocation  scattering  is  relatively 
higher. 

Finally,  Monte  Carlo  results  for  the  electron  mobility  as  a  function  of  the  2DEG 
density  are  shown  in  Fig.  30  for  operating  temperatures  of  300K  and  77  K.  A  monotonic 
decrease  in  mobility  with  channel  density  is  evident  at  both  temperatures  for  the  density 
range  chosen.  This  trend  is  in  keeping  with  experimental  observations*”’”^  and  results 
from  increased  carrier  scattering.  Electric  fields  at  the  hetero-interface  increase  with 
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Fig.  30.  Electron  mobility  vs  2DEG  density  for  temperatures  of  77  K  and  300  K. 

density,  and  hence,  carriers  are  more  strongly  confined.  This  leads  to  enhancements  in 
interface  roughness  scattering.  Furthermore,  the  interaction  with  optical  phonons  also 
increases,  but  is  not  as  strongly  dependent  on  density.  The  mobility  decrease  with 
density  is  the  largest  at  the  lower  77  K  temperature,  since  the  relative  impact  of  IRS  is 
stronger  than  at  300  K.  The  very  high  density  regime  has  not  been  probed  here.  It  is 
conceivable  that  at  high  densities,  the  Fermi  level  (Ep)  can  begin  to  exceed  the  sub-band 
energy  Ej.  This  could  then  change  the  phonon  scattering  magnitudes  as  Ep  -  E^  ~  ST© 
due  to  contributions  arising  from  the  phonon  emission  process.^^^  However,  the  behavior 
at  77  K  and  lower  temperatures  is  not  expected  to  be  affected,  and  the  trend  predicted 
here  should  remain. 

7.  SUMMARIZING  CONCLUSIONS 

Calculations  of  the  electronic  mobility  and  drift  velocity  have  been  carried  out  for 
bulk  GaN  and  AlGaN-GaN  heterojunctions  based  on  a  Monte  Carlo  approach.  The  bulk 
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calculations  were  intended  to  serve  as  a  validity  check  of  the  simulation  model.  For  the 
heterojunction  electron  mobility  calculations,  polarization  effects,  degeneracy,  and 
interface  roughness  scattering  were  all  taken  into  account.  Degeneracy  was  shown  to 
play  an  important  role,  especially  at  large  gate  bias.  Very  good  agreement  with  available 
experiments  was  obtained,  and  the  calculations  yield  a  set  of  best-fit  transport  parameters. 

Monte  Carlo  were  then  extended  for  calculations  of  the  large-signal  nonlinear 
response  characteristics  of  GaN-AlxGai.xN  High  Electron  Mobility  Transistors  (HEMTs) 
with  particular  emphasis  on  intermodulation  distortion  (IMD).  The  nonlinear  electrical 
transport  was  treated  on  first  principles,  all  scattering  mechanisms  included,  and  both 
memory  and  distributed  effects  built  into  the  model.  The  results  demonstrated  an  optimal 
operating  point  for  low  intermodulation  distortion  (IMD)  due  to  a  minima  in  the  IMD 
curve.  High-temperature  predictions  of  the  IMD  have  also  been  made,  by  carrying  out 
the  simulations  at  600  K.  An  increase  in  dynamic  range  with  temperature  has  been 
predicted,  due  to  a  relative  suppression  of  interface  roughness  scattering. 

Also,  a  model  analysis  of  the  large-signal  characteristics  of  GaN-AlxGai-xN  HEMTs 
was  carried  out  to  identify  parameters  and  their  optimal  settings  for  enhanced  microwave 
performance.  A  case  has  been  made  for  increased  mole  fraction  for  the  barrier  layer, 
reducing  the  transit  length,  and  introducing  a  thin  AIN  interfacial  layer  for  suppressing 
real  space  transfer  for  enhancing  the  device  performance. 

The  large-signal  response  characteristics  of  GaN-based  Gunn  oscillators  have  been 
re-examined  based  on  recent  bandstructure  data  that  has  become  available  for  the 
wurtzite  material.  Wurtzire  GaN  is  also  more  interesting  since  it  has  been  shown  to  have 
a  higher  drift  velocity  that  should  promote  higher  frequencies  for  Gunn  diode  operation. 
Device  performance  parameters  of  interest  were  calculated  for  two  separate  situations. 
First,  devices  with  the  conventional  single  notch  structure  were  simulated  to  obtain 
predictions  of  the  output  power  and  conversion  efficiency.  Then,  multiple  domain 
devices  containing  repetitive  serial  segments,  were  analyzed.  The  repetitive  structure  has 
been  studied  here  since,  to  the  best  of  our  knowledge,  the  benefits  and  potential  of 
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multiple  domain  structures  (MDS)  have  not  been  examined  in  the  context  of  GaN  Gunn 
diodes.  Theoretical  predictions  for  GaAs  Gunn  oscillators  call  for  an  scaling  in  output 
power  with  the  number  of  segments  N. 

Device  parameters  such  as  the  operating  frequency,  power  output  and  conversion 
efficiency  were  obtained.  Variations  due  to  changes  in  temperature,  biasing  voltage  and 
device  length  were  also  included.  As  expected,  the  operating  frequency  has  been  shown 
to  fall  with  increasing  temperature  due  to  a  lower  drift  velocity.  The  results,  for  a  single 
transit  region  structure,  showed  that  conversion  efficiencies  up  to  about  2  percent  can  be 
attained  at  300K  and  would  be  around  1.45  percent  at  450  Kelvin.  The  room  temperature 
value  is  higher  than  that  reports  for  GaAs  devices.  Also,  the  efficiencies  are  slightly 
higher  than  those  calculated  elsewhere  for  GaN  on  the  basis  of  hydrodynamic  models. 
Corresponding  peak  power  outputs  of  0.68  W  and  0.52  W  have  been  predicted  at  300  K 
and  450  K,  respectively.  Here  again,  considerable  improvement  is  predicted  over  other 
material  (e.g.  the  reported  results  of  about  0.150  W  for  InP.)  However,  it  must  be 
emphasized  that  since  device  self-heating  effects  have  not  been  included  here,  these 
values  represent  an  overestimation. 

For  multiple-domain  structures,  considerable  improvements  in  output  power  have 
been  predicted.  Without  changing  the  device  area,  the  output  power  results  were  seen  to 
roughly  scale  with  the  number  of  segments  N.  Slight  deviations  (sub-linear  behavior) 
seen  can  be  attributed  to  the  loss  in  absolute  repetitiveness  from  one  segment  to  the  next. 
Hot  electron  effects  and  larger  average  carrier  energies  contribute  to  the  sub-linear 
behavior.  Hence,  for  devices  with  with  scaled  increases  in  area,  the  output  power  would 
not  be  quite  as  large  as  the  predicted  factor.  Also,  the  noise  for  the  multi-domain 
devices  is  expected  to  be  higher  than  the  conventional  structures.  This  factor  has  not 
been  studied  and  will  be  analyzed  elsewhere.  In  any  case,  the  advantage  of  multiple  GaN 
Gunn  diode  structures  has  been  made  obvious,  and  such  device  do  merit  serious 
experimental  study. 
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Monte  Carlo  studies  of  2DEG  transport  were  used  to  analyze  carrier  mobilities  in 
wurtzite  GaN  HEMTs  taking  account  of  the  interface  polarization  charges  arising  from 
edge  dislocations  strains  and  all  other  pertinent  carrier-phonon  scattering  processes. 
Dislocation  scattering  has  been  shown  to  be  important  in  bulk  GaN,  but  has  not  been 
extensively  studied  in  the  context  of  GaN  HEMTs.  Explicit  expressions  for  dislocation 
scattering  were  developed  in  the  context  of  weak  perturbation  theory  in  the  Born 
approximation.  It  was  shown  that  dislocation  scattering  is  anisotropic,  and  would  have  a 
much  stronger  influence  at  lower  temperatures.  For  self-consistency,  details  of  the  carrier 
density,  sub-band  energies  and  eigenfunctions  were  included.  This  was  done  by 
including  two  equivalent  1x1  Hamiltonians  for  the  37c  conduction  bands,  and  a  6  x  6 
Hamiltonian  for  the  A,  B,  and  C  valence  bands.  The  eigen-value  information  was  used  to 
calculate  the  scattering  matrix  elements  using  parameterized  envelope  wavefunctions. 

Finally,  electron  mobilities  were  calculated  as  a  function  of  operating  temperature  for 
various  dislocation  densities.  The  trend  was  shown  to  be  in  keeping  with  the 
experimental  reports.  Room  temperature  electron  mobilities  were  predicted  to  be  around 
1.711  X  10^  cm^A^s,  which  is  fairly  close  agreement  with  the  data  reported  by  the  Santa 
Barbara  group^°’^^^.  At  the  highest  dislocation  density  of  10^°  cm“^  ,  mobility  reduction 
of  16.8  percent  and  8.6  percent  are  predicted  for  77  K  and  300  K,  respectively.  Given  the 
sharp  drop  in  phonon  scattering  at  the  lower  temperature,  the  impact  of  dislocation 
scattering  is  relatively  higher.  Monte  Carlo  results  for  the  electron  mobility  as  a  function 
of  the  2DEG  density  were  also  obtained  for  operating  temperatures  of  300K  and  77  K. 
The  mobility  decreased  monotonically  due  to  increased  interface  roughness  scattering. 
This  trend  is  again  in  keeping  with  experimental  observations®®’^^^.  The  mobility  decrease 
with  density  was  the  largest  at  the  lower  77  K  temperature,  since  the  relative  role  of  IRS 
is  stronger  than  at  300  K. 
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